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ABSTRACT 


This thesis presents the development and validation of 
a traffic circle simulatcr for use as a tool in the 
evaluation of traffic circle performance. The simulator was 
constructed as a general traffic flow model. Different 
traffic circle situations can be simulated by supplying the 
Simulator with the descriptive geometrical characteristics 
of the circle in terms cf circle configuration, and 
pertinent data including the total traffic volumes from the 
various approaches of the circle. 

The model has the capacity for the dynamic analysis of 
traffic flow at a circle with up to 6 arteries and up to 2 
approach lanes, 2 exit lanes and 2 circle lanes. The model 
was implemented on an IBM 360/767 computer in FORTRAN IV. 

The model was validated at the macroscopic level. 
Several different macroscopic comparisons were made between 
simulated phencmena and real data collected at various 
traffic circles in the City of Edmonton. The comparisons 


were found to be consistent and reasonable. 
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Intreduction 


1.1 General Introduction 

The increased use of motor vehicles for the past three 
decades, and the prospects cf even greater increases in the 
future, have alerted traffic engineers in particular and 
other public cfficials in general to the problem of 
providing adequate facilities for safe and efficient 
movement of traffic on the roads. The problems in the 
cities and metropolitan areas have already reached major 
proportions. Freeway develcpments, which emerged as a_ sure 
answer to the urban and metropolitan traffic problems, are 
becoming less attractive as a means of accomodating travel, 
due to extremely high cost and adverse social effects. 
Therefore more emphasis is now being placed on obtaining 
maximum capacity and efficiency of the existing street 
networks as well as improved public transit. iieesuch eed 
direction, detailed study should be given each element, to 
gain complete understanding of its operation and its 
relation to each of the other elements in the system. 

Every driver would like to proceed as he pleases 
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through the street network from his origin to his 
destination. Since his path crosses that of other vehicles 
at intersection points in the system, it is desirable and in 
fact imperative to minimize the chances that the 
intersection of the vehicle paths will not result in 
collision or undue delay to any one vehicle or group of 
vehicles. 

The paths of two vehicles can be separated by either 
Space or time. When a few vehicles are distributed over a 
relatively large number cf streets, as happens in rural 
networks, the separation in time due to the low probability 
of interference usually cbviates any need to control the 
intersection points. When street use becomes more intense, 
however, the probability of time separation of vehicles by 
mere chance beccmes smaller, and the possibility of vehicle 
collisions at potential intersections becomes a concern to 
public officials. When the severity of the problem 
justifies ist, vehicles can either be separated at 
intersection points by space (grade Separation) or rules can 
be established to force time separation. 

The rules applied vary from the simple right-of-way 
rule to traffic control devices which range from ‘yield 
signs' to ‘traffic signals'. The degree of control should 
increase with increasing probability that time separation of 
the vehicles at intersecticn points will not occur by mere 


chance. Therefore the rules change from those requiring 
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restriction of cne or two streams ( YIELD and TWO-WAY-STOP) 


to those requiring restriction cf all intersecting streams 


( FOUR-WAY-STOP, TRAFFIC CIRCLES, SIGNAL CONTROL). 
question to ask at this juncture, is how to organize 
control traffic so that drivers could fully enjcy 
increased mobility on the reads. 
"To answer this question is, in essence the job of 
traffic engineer. Ccnsidering such objectives as 


minimization of delay or travel time, a decrease 
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the number of stops, an increase in the smoothness of 


driving, and an increase in safety, he must 


traffic to the existing road system by regulation 


fat 


and 


control or he must fit the road system to the traffic 


demands by planning and redesigning."{ 14; p. 2]. 


1.2 Nature of the Problen 


The flow cf traffic through the urban and metropolitan 


areas could be thought of as vehicles moving on roadways 
intersections cf readways under the following methods 
devices of intersection ccntrol: 

a) Yield signs; 


b) Stop signs; 


c) Traffic circles; 
d) Fixed-time traffic signals; 
e) Vehicle-actuated signals; and 


f£) Computer-controlled signal-systems. 


and 


and 
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Until a vehicle encounters an intersection (where it might 
experience some delay), it is assumed that the vehicle 
travels at a desired speed unless it follows a_ slower 
vehicle. Again, a vehicle could be delayed when it is 
forced to travel below the desired velocity due to a traffic 
conjestion downstream occassionally caused by a restriction 
SEA hS roadway such as construction, accident, narrowing and 
sO on. We could say therefore that vehicular delays are 
caused mainly by the vehicles encountering intersections 
controlled by any of the six methods mentioned above. 

While delay is undesirable from the driver's point of 
view, it serves some form of useful purpose from the point 
of view of the road designer or the traffic engineer. Delay 
is useful in describing the level of service at an 
intersection or ina system of streets. Also it lends 
itself to economic analysis. The driver, however, is 
annoyed by delay, and he is constantly attempting to 
eliminate it altogether cr at least minimize it. It is 
useful therefore, to know what the major causes of delay 
are, andto what extent each contributes to the total delay 
in a complete network of streets. With such a knowledge, 
the engineer can evaluate what the overall effect would be 
with the removal of one of the causes of delay. 

With the above reasoning, it is important therefore 
that the traffic engineer knows which method of control is 


best for a given intersecticn condition, having taken into 
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account the volume of traffic at the intersection, and the 
turning volumes from the various approach streams of the 
intersections. However, only meager informaticn is 
available concerning methcds of control below the level of 


traffic signals. 


1.3 The Need for Traffic Circle Simulation Model 

The overall cbjective of this thesis is to better 
identify the operating characteristics that affect the level 
of service at a traffic circle, and in this way give the 
traffic engineer some criteria for the measure of 
performance at a traffic circle. Once this type of 
knowledge has been gained, it will be relatively easy to 
develop criteria for deciding which method of control to 
apply to a specific intersection (stop signs are usually not 
of concern, since they are converted to signals when 
accidents or traffic vclumes warrant such installations). 

Most often the traffic engineer has available to hin 
some form of traffic data obtained from field surveys. This 
data is usually made up of the total volume inflow of 
vehicles from the various approach directions of the 
intersection and the turning volumes from same (see Table 
1.1). Table 1.1 is a summary of morning 15 minute peak 
period traffic counts stratified into the various turning 


volumes for cars and trucks separately. 
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Briefly, the abbreviations in the Table mean the following: 


280 ST = Code name fcr Groat Road; 
C = Cars; 

aT = Trucks; 

CRSSWLK = Crosswalk; 

PEDST = Pedestrian; 


NTH, STH, EST, and WST respectively stand for North, 

South, East and West. 

Presently, there are many simulation models for signal 
controlled intersections, so that the engineer could always 
apply his data to those mcdels; modify the data as he 
wished, to study the intersection performance under various 
operating characteristics [13, 17, 29, 31]. It cannot be 
overemphasized therefore, that a similar simulation model on 
the traffic circle should be available to the traffic 
engineer since the flow of vehicles through a _ circle 
deviates very much from the passage of vehicles through an 
orthogonal intersection under any other form of ccntrol. 

During the passage of vehicles through a _ traffic 
circle controlled by unsignalized devices such as the YIELD 
signs, vehicles cannct be assigned a generalized velocity 
such as the "constant queue discharge velocity’ (a 
hypothetical velocity which describes the departure of 
queued vehicles at a signalized intersection during the 
green phase). This is because each individual vehicle which 


queues up at the entry to the circle has to evaluate and 
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accept or reject gaps in the circle traffic stream before 
entering the circle. The passage of vehicles through a 
traffic circle could therefore be likened to the passage of 
a minor street vehicle through its intersection with a major 
street. But, as compared to the traffic circle, the passage 
of such a minor street vehicle through the intersection 
could be thought of as being instantaneous (in the sense 
that, once it accepts a lag or _ gap to enter the 
intersection, alee does not spend much time at _ the 
intersection). 

Operational rules at a traffic circle require more 
interpretation on the part cf the driver, and more activity, 
Skill and alertness than at a signalized intersection. In 
order to study the control of traffic behaviour at traffic 
circles, an attempt must be made to consider fully the 
behaviour aspect of the operation, since there is more 
interplay between drivers at traffic circles than at 
Signalized intersections. Another important characteristic 
of a traffic circle is the fact that the rate and time of 
arrivals of vehicles is of major importance, because it is 
vehicle presence in the circle stream that determines what 
the driver in the approach stream does. With the traffic 
circle simulation model the traffic engineer could study the 
performance cf the circle through the use of his field data 
Or thrcugh arbitrary apportionment of generating 


characteristics to the various approaches of the circle. 
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Presently, there is a great concern in urban areas in 
most countries especially in North America, where the use of 
traffic circle is very limited, to replace existing traffic 
circles with signal controlled intersections and to halt the 
construction of more traffic circles. But such statements 
are usually general policy statements, since there is 
normally limited valid evidence to Warrant such 
undertakings. 

Before a traffic engineer decides to replace an 
existing traffic circle with a Signal-controlled 
intersection or vice versa, he must have subjected the 
particular intersection to variable traffic conditions 
through the use of a signal-controlled intersection model 
and a traffic circle simulaticn model. Then, he must have 
convinced himself that improvements of the intersection 
performance expected through such a conversion would be 
worth the efforts involved in the implementation of such a 
decision, especially before committing public funds. 

Of course, such an evaluation process can only be 
carried out cheaply, realistically and conveniently through 
simulation. He needs Simulation models capable of 
simulating the operations of the intersection under signal 
control and under traffic circle control. He might arrive 
at the conclusicn that the volume of left-turning traffic 
from a particular approach would make a traffic circle 


superior to a signal controlled intersection even with 
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*protected* left turns, or vice versa. In short, with the 
capability of varying input information to both the 
Simulation models, the engineer could evaluate the 
justification of any type of conversion from the use of one 
method of contrcl to the other. 

The above analysis of the problem of implementing some 
conversion policies, greatly justifies the need for the 
traffic circle simulation model as described in this thesis; 
a model which accepts traffic data of all forms, dynamically 
varies the arrival frequencies and turning movements during 
a stipulated period of the simulation, at rates specified by 
the traffic engineer. A model with such a capability will 
therefore give a good indication of the real life 
performance expected from the traffic circle under the 
specified conditions. 

This thesis presents the development and validation of 
such a traffic circle simulation model for use as a tool in 
the evaluation of traffic circle performance in terms of 
level of service. Chapter II is a general review of 
existing literature covering the major analytical and 
Simulaticn work done in the field of Traffic Flow. Chapter 
III describes "TRACISM' the Traffic Circle Simulation Model. 
The Chapter is devoted to the formulation of the traffic 
circle simulaticn model. The important parts of the model 
are discussed with the relevant diagrams. The detailed 


analysis of the representaticns of the model in the computer 
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and the actual flow logic are described in Appendices E and 
C respectively. 

Chapter IV describes the validation studies cf the 
model by a series of simulation runs on test nodel 
configurations, and the eventual validation of the model by 
comparison with real data. A major part of this Chapter 
deals with the study of various measures of effectiveness of 
the circle performance. Chapter V contains the various 
conclusions as made evident from Chapters III and IV and 
recommendations for further study on the Traffic Circle and 


their role in urban street networks. 
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Chapter Ii 


Traffic Flow Theory and Simulation 


2.1 Introduction 

With the recent concern for the control of urban air 
pollution, vehicle manufacturers are urged to produce 
"roadworthy" cars with respect to pollution. In much the 
same wWay, automobile users are concerned about the roads on 
which they drive, and therefore urge the responsible public 
officials to provide ‘carworthy* roads. Carworthy roads 
need not only be smooth with respect to potholes but also 
should afford smooth trifs. This can be achieved through 
redesigning and control of the existing roads. 

The importance attached to the control of traffic on 
the roads is evident from the numerous attempts that have 
been made, through the developements of various traffic 
models aimed at improving cne traffic situation cr ancther. 
Insights into the behaviour of traffic on the roads are 
obtained through two major approaches, namely: 

a) Experiments and field observations and 

b) Modeling of the real situation. 

The former methcd is not entirely satisfactory. Conclusions 
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drawn from field cbservations may be wrong because real 
traffic conditions usually make it impossible to reproduce 
the causes of the observed behaviour. The field experiments 
are extremely useful however, and should be conducted before 
a large capital outlay is made since they give an indication 
of the validity of the latter method. Since modeling is 
essentially experimenting on a working analogy having the 
necessary relations and properties similar to those of the 
system being studied, modeling techniques may be classified 
as follows: 
(i) Analytical model. 
In this case equations sufficiently representative of 
the system are solved analytically. 
(ii) Software simulation model. 
Here a large digital computer may be used to assist in 
the solution. This approach is particularly relevant 
if nonlinearities are present in the system. 
(iii) Hardware simulation mcdel. 
Here a special-purpose computer constructed from units 
performing mathematical and/or logical operaticns on 
simulated variables is employed. 
The study of the system behaviour using the software 
Simulaticn model or the hardware simulation model is 
hereafter referred to as ‘Simulation’. The next three 
sections give a summary of analytical models in traffic flow 


and the simulaticn cf traffic flow. 
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2.2 Analytical Models in Traffic Flow 

Solutions to traffic flew problems have been attempted 
by analytical methods for many _ years. The analytical 
approach has’ keen confined mainly to analysis of 
uncontrolled intersections, intersections controlled by 
yield signs or two-way stops, intersections controlled by 
fixed-cycle traffic signal and those controlled by vehicle- 
actuated traffic signals. Most of the analytical models on 
the fixed-cycle traffic signals, have been concerned with 
the queue measurements at the beginning of red periods. The 
fixed-cycle traffic signal queues have been investigated by 
probabilistic methods by a number of investigators, anda 
brief review of the various methods of approach is given 
below. 

Beckman, McGuire and Winsten [3] considered a discrete 
time queuing model with binomial arrivals and regular 
departure headways during the green phase. They went cn to 
derive a relation between the stationary mean delay per 
vehicle and the stationary mean queue-length per vehicle at 
the beginning of a red period of the traffic signal. Haight 
[24] and Buckley and Wheeler [8] considered models with 
Pcisson arrivals and again regular departure headways during 
green phase, and investigated certain properties of the 
queue length. Newell {42] dealt with the model proposed by 
Beckman et al, and obtained the probability generating 


function of the stationary queue length distribution. 
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Darroch {11} discussed a more general discrete time model 
with stationary, independent arrivals and regular departure 
headways, and derived a necessary and sufficient condition 
for the staticnary queue length to exist, and obtained its 
probability generating function. 

McNeil [36] proposed a model for the solution to the 
fixed-cycle traffic signal problem for compound Poisson 
arrivals. In McNeil's model, an intersection is controlled 
by a traffic signal with fixed cycle time T, the possibility 
of other delays such as those due to turning vehicles, is 
completely ignored. Arrivals at the signal form a compound 
Peisson process, which, as McNeil himself points out, is 
unrealistic since such an assumed arrival process presumes a 
zero car length. If vehicles arrive to a green signal and 
an empty queue they cross the intersection undelayed, while 
on the contrary, the vehicles depart when they reach the 
head of the queue, provided the signal is green, each 
vehicle taking a constant time to move off. McNeil's model 
contains numerous assumptions typically embodied in most 
analytical models aimed at simplifying the situation under 
study enough to make it mathematically tractable. The 
assumptions usually render the models too simple and 
unrealistic to be useful. 

Vehicle-actuated signal intersections have been 
investigated in two ways, the non-priority and the priority 


types. In the non-priority type, neither of the 
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intersecting streets is favcured, rather the intersection is 
controlled by signal which favours one queue at a time and 
switches when the favoured queue empties. The priority 
vehicle-actuated traffic signal on the other hand is the one 
at which the main street signal is green unless traffic 
arrives on the side-street. 

In the priority vehicle-actuated traffic problem, 
Little [33] has considered a queuing process of side street 
traffic. The intersection considered is made up of two one- 
way streets. But there is no apparent reason why the 
intersecting streets should be assumed to be one way other 
than to make the mathematical analysis simpler, and of 
course to limit the scope of the problem. In Little's 
model, the signal shows green on the main (priority) street 
unless it is actuated to change by side-street traffic. 
There is no detector on the main street so that the density 
of the main-street traffic has no effect on the operation of 
the signal. 

Side-street traffic causes the signal to change by 
means of a detector which senses whether or not there is a 
vehicle at the intersection. After detection of a vehicle, 
an activation period t, must elapse before the Signal 
changes to green for the side-street traffic, this period 
corresponding to the amber (clearance) period cn the main 
road. After changing, the signal remains green to _ side- 


street traffic for a fixed interval of length g, and each 
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green period is followed by a red period of minimum length 
cr. During this minimum red period, the length of which isa 
constant, programmed into the signal mechanism, the signal 
will not change tc green for the side-street traffic. 
Little's model is very similar to the queue balancing model 
of Lieberman [32], where the signals dc not respond to 
vehicle presence but rather to some threshold values of 
queue lengths and signal lengths. Little develored a 
mathematical model of the intersection operation as 
described above and obtained the following: 
1. The transition probabilities of the queue lengths at 
the end of each cycle (green plus red); 
2. The probability distribution of the queue lengths in 
the steady state; and 


3. The mean gueuve lengths. 


2.3 Simulation 

Simulation could be defined as a dynamic 
representation of reality, achieved by building a model and 
moving it threugh time. The problems solved by simulation 
are characterized by being mathematically intractable and 
having resisted solutions by analytical methods. The 
problems usually involve many variables, many parameters, 
functions which are not well-behaved mathematically, and 
random variables. Simulation should therefore be the 


technique to apply when the number of nonlinearities in a 
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system makes analytical techniques inapplicable. After 
identifying and formulating the problem, the next major step 
in simulation is the modeling phase which is described 


below. 


2.3.1 Modeling and Computer Simulation 

The follcwing concert of modeling and simulation is 
taken from C.W. Bell and R.N. Linebarger [4] who define 
modeling and simulation as consisting of four sequential 
steps as follows: 

a) Concept development; 

b) Model formulation; 

Cc) Model implementation; 

d) Model operation. 

The ability to simulate system behaviour depends upon the 
formulation cf the model, and each step in the modeling and 
Simulaticn process contains a feedback path whereby 
implementation and evaluaticn of current models can generate 
new structures and question of validity. 

The first step involves the development of the 
conceptual system behaviour. Here, the analyst observes and 
classifies phenomena into ccnceptual framework. The second 
step is the model formulation which may be done in the 
following four ways: 

(i) Verbal; 


(1i) Physical; 
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(iii) Mathematical. 

The verbal model provides the structuring of the system 
concepts on a verbal basis and seldom undergoes 
implementation cf any kind. Rather, the verbal model is 
primarily used for developing new concepts of system 
behaviour. The physical model is usually a scaled analogy 
of the actual system under study. 

In mathematical models the system to be studied is 
represented by mathematical relations derived from 
observations of fundamental data. The mathematical model is 
almost always manipulated by means of computer simulation of 
the basic system state-transition functions. 

After formulation, the next process is implementation. 


Here the abstract model is translated into some fcrm of 


device which allows the mcdel parameters and system 
structures to be manipulated. Physical models are 
implemented using physical analogies. Mathematical models 


are implemented using a variety of tools, primarily 
computer-oriented. 

The final step in the simulation process is the 
operation of the implemented model. Here, the model 
parameters can be altered or varied and the resultant model 
response analysed, On the basis of these results, 
refinements can be proposed for both the formulated and the 
implemented model, to make them better representations of 


the system being studied. 
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In summary, the process of modeling and simulation 
involves a set of sequential steps with feedback paths. 
Inherent in computer simulation studies is the principle of 
interaction with an implemented model whereby the model can 
be exercised and the resultant response compared with 


observed data. 


2.3.2 Analytical versus Simulation Models 

There are several Similarities and differences between 
analytical models and simulation models. Scme of the 
Similarities are: 

1. Both appreaches require a thorough understanding of 
the process being mcdelled. 

2. Both appreaches involve the formulation of an abstract 
model which represents a concrete situation. 

3. The translation of a real situation into either type 
of model reflects the users concept of what the key 
elements cf the system are and how they interact. 

Some of the major differences between analytical and 
Simulaticn models are: 

1. The analytical model yields analytical solutions 
whereas the simulaticn model shows what happens under 
a particular set of assumptions and does not yield a 
*solution*’. 

Py An analytical model, in order to be solvable is often 


so gross a simplification of the actual situation as 
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to yield invalid results. Simulation usually permits 
a less abstract and relatively more faithful 


representaticn of a real system. 


2.3.3 Traffic Simulation Models 


As applied to traffic flow studies, the following is a 


summary of the requirements of a good simulation model as 


listed by Drew [12]: 


a) 


b) 


Cc) 


d) 


e) 


f) 


It must provide an easy, inexpensive method of traffic 
Simulation. 

It must ke general enough so that any Similar traffic 
configuration can be simulated by the input of the 
proper geometric parameters. 

The input to such a program must be easily understood 
and capable of execution by non-computer oriented 
personnel. 

It must furnish output which is easily readable and 
which contains all parameters needed by the traffic 
engineer for application in the design or modification 
of traffic streams. 

It must be written in modular fashion, such that any 
of the moduli can be changed without affecting the 
rest of the program (for example, the car-following 
process should be completely independent of the input 
generation process, ard so on). 


It must be written such that new moduli such as 
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traffic hazards, curves, and some other variable 
traffic conditions can be added without extensive 
programming changes. 

It must be machine independent, written in one of the 
higher level languages such as FORTRAN IV in such a 
Manner that a relatively unsophisticated programmer 


can modify it. 


2.3.4 Advantages and disadvantages of Traffic Simulation. 


a) Advantages: As applied to traffic flow studies, 


the following is a summary of the advantages of simulation 


again listed by Drew [12]. 


‘be 


The task cf laying out and operating a simulation is a 
good way of systematically gathering pertinent data. 
It makes for a broad education in traffic 
characteristics and oferation. 

Simulation of complex traffic operation may provide an 
indication of which variables are important and how 
they relate. This may lead eventually to successful 
analytic formulaticns. 

Simulation gives an intuitive feeling for the traffic 
system being studied, and is therefore instructive. 
Simulation gives a control over time. 

Real time can be compressed, and the results of a long 
signal phase can be cbserved in a few minutes of 


computer time. On the other hand, machine time can be 


aolisidnte ot22ear 30 eegniasvbdalh bap eBGepaevbs Pal sh 
,sebbote wot? oittes? of betiqgs 2&4 rzapsrmevbs = (s 
soldstieie io eepstdevbs. si? to veneers B 22 patwollco? od? 
[Sta wort ye a 
s et ooitsiueie 5 paitateqo fas tuo patyed to. #e59- oT - 
cjed stonteiag gabxedrey ylisaltemedaya to yaw Soop 
oitaert al adttsowhs Bsoxd) @ 302 seeded oF 
.aoltszego bas exivetzetoombio 
ne shivosy (ee voltsaeqo Ti32sat xelqaon 20 conan a | 
wod to5 stnasttoqeit sts asidetisy dotty 20. apedeokbak | 
ipestsoous of pilsdonors Brot yaw 2idt setstes ved? 
ofltars oft re notretaake «6 


b) 


23 


expanded and run slower than real time so that all the 

manifestations of the complex interactions of traffic 

movements can be comprehended. 

Simulation is safe. It provides a means of studying 

effect of contrcl measures on existing road systems. 
Disadvantages: 


There are however some disadvantages to the 


application of simulation to traffic flow studies. The 


major disadvantages are listed by Ferris [14] as follows: 


1. 


Considerable effort is required in the programming of 
a Simulaticn model of some desired complexity. 
Sinulation of traffic flow problems becomes’ very 
costly, and consideration must be given to man-machine 
feasibility. 

Accurate and detailed data is essential to any traffic 


Simulation model. 


However, it should be mentioned that some of the above 


disadvantages apply equally to analytical models, so that 


the 


utility of simulation models is still superior to that 


of analytical mcdels. 


2.3.5 Languages used in Simulation 


The languages used in Computer simulation studies fall 


into two major categories, namely: 


a) General Purpose Simulation Languages, and 


b) General Purpcse Programming Languages. 
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The early 1960's saw the advent of several general purpose 
Simulation languages. Most of them were adaptaticns of 
FORTRAN (for example, SIMSCRIPT, DYNAMO, GASP) or ALGOL (for 


example, SIMULA). Some were entirely new, for example GPSS. 


a) Advantages: 

1. The main advantage of the general purpose simulation 
languages as compared to the ordinary general purpose 
programming languages is the fact that the fcrmer 
incorporate means fcr controlling the sequence in 
which events occur. This sequencing aspect of 
Simulation models introduces many complexities when 
Ordinary programming languages are used, in which case 
the co-ordination of events in the model remains the 
responsibility of the analyst or sometimes the user. 

2. Another major advantage of simulation languages is 
that their diagnostic programs can check for logical 
errors as well as syntax and capacity violations in 
the model. 

3. Simulation languages usually have built-in provision 
for collecting and printing out at least some of the 
statistical outputs desired from most Simulation 
models. 

4. Finally the use of simulation languages makes the 
model closer in outward appearance and structure to 
the situation being simulated than when nonsimulation 


language are used. 
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b) Disadvantages: 
1. The major disadvantage of most general purpose 
Simulation languages is the fact that they require a 
large-scale computer and are often therefore limited 
in use. For example, GPSS/360 and Flow Simulator are 
designed for computers with core capacity of 64K or 
greater. 
2. Another disadvantage of the use of most general 
purpose Simulation languages is the fact that they are 
often too slow for large scale simulation systems. 
Presently, there are dozens of simulation languages, 
but most of them are special-purpose types which are not 
widely used. The major simulation languages in North 
America today include GPSS, SIMSCRIPT, SIMULA, and DYNAMO. 
Of these, GPSS is the most widely used. GPSS has been more 
widely employed in traffic flow simulation studies than any 
of the other simulation languages. Before the advent cf the 
general purpose simulation languages, general purpose 
programming languages such as FORTRAN or ALGOL had been used 
for many years for simulation studies and are still in use. 
Presently, the most widely used nonsimulation languages are 
FORTRAN, ALGOL and PL/1. Of these FORTRAN is used more 
often than the rest. 

Watjen {51] used GPSS for the simulation of a_ simple 
traffic network, tc study the effect of offsets in signal 


light synchronization. A pioneer application of GPSS to 
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the simulaticn of traffic flow was done by Blum [6]. He 
used GPSS to simulate a network of intersections where the 
intersection area is divided up into cells represented by 
facilities in GFSS. 

A more recent application of GPSS to traffic flow 
Simulaticn studies has been done by Lieberman [32] at the 
General Applied Sciences Lakoratories,Inc. in New York. The 
model is an integrated set of computer programs which 
Simulate traffic flow at signal-controlled intersections. 
The complete set of programs are collectively named SURF 

(Simulation of URban traffic Flow). 

The intersection area in the model is subdivided into 
facilities some of which are accessible to both pedestrian 
and vehicular traffic. Separate queues are alloted for 
vehicular and pedestrian traffic. Each vehicle is treated 
as a separate transaction whose behaviour is monitored from 
the instant it is generated into the system until it crosses 
a terminate line and is discharged from the systen. 

The entire approach distance of each approach lane 
represents a gueue, whereas the distance from a stop line to 
an intersection entry line of each approach lane represents 
a facility which is used by both vehicular and pedestrian 
traffic. Similarly, each distance of an exit lane from the 
intersection exit line to the terminate line (concurrent 
with the stop line at an approach) constitutes a facility 


which is again used by both vehicular and pedestrian 
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traffic. Thus, upon entry into the system, a vehicle 
travels the approach distance by spending some calculated 
time in the appropriate gueue. Upon emergence frem the 
approach queue, the vehicle begins its journey into and out 
of the intersection by successfully seizing (occupying) and 
releasing (leaving) facilities, until 2c has finally 
released a terminal facility after which it leaves the 
systen. 
Some prime features that are exhibited by Lieberman's 
microscopic model are: 
1. The ability to accomodate a wide variety of 
intersecticn configurations; 
2. The ability to accept multiphased traffic signals; 
3. The ability to simulate traffic-actuated signals; and 
4. The ability to accomodate coupled pedestrian and 
vehicular traffic. 
Vehicle arrival is generated by approach according to 
a shifted exponential distribution with a minimum inter- 
arrival gap of one second, and a mean gap equal to the 
reciprocal of the specified flow rate (traffic volume in 
vehicles per hour). Turning maneuver is assigned 
stochastically based upon the specified turning percentages 
for each approach. A turning vehicle joins the appropriate 
turning lane at the approach and a through vehicle joins the 
shorter of the approach queues. 


This methcd of vehicle arrival generation and eventual 
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assignment to a lane depending on turning maneuver of a 
vehicle, eliminates the possibility for lane-changing in the 
approach flow. This is because, during peak hours, lane- 
changing within the vicinity of the intersection is almost 
always impossible, so that upon entry into the approach, 
drivers quickly determine their lanes and stay in that lane 
for intersection maneuver. 

A major objective of the simulation was to study the 
difference in intersection performance under fixed-cycle 
Signal control and varying cycle length signal control. To 
this end, the following test models were used for validation 
purposes. 

Model A: Two approach lanes in each direction with 
traffic controlled by a two-phased signal of fixed cycle 
length. 

Model Bs: Iwo approach lanes in each direction with 
traffic controlled by a two-phased signal of varying cycle 
length. The signal timing is controlled by the density of 
traffic flow. 

Incorporated in Model B is a well-defined algorithm of 
queue-balancing strategy employed by Lieberman. The 
objectives and the consequent measures taken to achieve the 
objectives through the use cf the ‘responsive traffic 
control algorithm" are: 

1. Optimization of intersection performance by attempting 
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a) utilize intersection capacity fully; 
b) reduce total delay time (in queues) to a 
hinimun; 
c) equalize (in each direction) the queue 
lengths and delay times. 
2. Reduction of driver irritation and consequent accident 
potential with the following policies: 
a) Set maximum signal phase length; 
b) Set minimum Signal phase length. 
3. Prevention of spill-back into "upstrean' 
intersections. 

With the akove objective, the following specifications 

were made: minimum signal length, maximum signal 

length, amber signal length, critical queue length for 

East-West traffic (EW) and critical queue length for 

North-South (NS) traffic. 

During the simulation, traffic signal ‘interrogates! 
intersection traffic status every Simulation cycle (scan 
interval) and acts according to prevailing conditions. A 
signal remains green in a direction such as East-West unless 
any of the following conditions is met. 

a) Maximum signal for the direction is exceeded. 
b) If signal has been green longer than the minimum 

specified and not more than one vehicle occupies a 

queue in this direction. That is, if demand is 


dissipated in this direction. 
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c) If signal has keen green longer than minimum length 
Specified and critical queue length in the other 
direction has been attained or exceeded, while the 
queue length in this direction is shorter than the 
prescribed critical value. 


) meet thes > oie Evs ¢ where 


dc average length of North-South queues; 
Cie average length cf East-West queues; 
E..= net inflow cf vehicles in North-South lanes 


during period cf minimum signal length. 

Hence, the signal remains green until the queves in the 
other direction become slightly longer. For the Model A, 
fixed specifications were made for Green, Amber and Red for 
each direction. The statistics which were generated by the 
program included the following: 

(i) Variation of queue lengths as a function of time; 

(ii) Travel and delay times for vehicles according to 

discharge lane and turning maneuvers; and 

(iii) Travel and delay times for pedestrians as well. 
Other Ststastics included the intersection capacity, 
expressed in vehicles per hcur. After subjecting the two 
models to various rat foc conddsbionsig(ot, tie typical 
intersection configuration, Lieberman [32] reached the 
following conclusions: 

1.. There is no advantage from the use of a computer 


controlled signal during periods of saturated demand 
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when the demand conditions exceed the intersection 
Capacity. 
2. When demand conditions exceed intersection Capacity, 
there may be only slight improvement in overall delay 
time when a responsive signal replaces a fixed time 
Signal. But the responsive signal gives a remarkable 
improvement in the balancing of queues at the various 
approaches. 
3. Intersection performance improves subtantially in all 
respects, if the incoming flow rate is at, or below 
saturation level when a computerized signal system 
replaces a fixed cycle system. 
A novel feature in Lieberman's model is its ability to 
simulate varying (peak-hcur) demand conditions. When 
varying traffic demands are simulated, queue lengths and 
delay times are greatly reduced during peak hour traffic 
conditions due to the response of the signal to changing 
traffic demand in the use of the responsive systen. 

Lieberman's model is very realistic and excellent as a 
tool for helping the traffic engineer in the evaluation of 
intersection performance. The model is more of a tool than 
most other models because in Lieberman's model the traffic 
engineer does not have to worry about input data such as 
vehicle acceleration, vehicle deceleration and other such 
"hard-to-collect' data. 


The input data required by the system is precisely 
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that information which is available to the traffic engineer. 
This data consist cf: 
a) density of vehicular and pedestrian traffic flow from 
each direction; 
b) vehicular turning mcvements from each direction; 
Cc) signal phasing and cycle time (or actuated control 
specifications). 
This description of Lieberman's simulation program shows how 


GPSS can be applied in complex traffic situations. 


2.4 Simulation of Traffic Flow 


2.4.1 Introduction 

To fully evaluate the effect of a traffic situation, 
it would be desirable to test it under a wide variety of 
conditions. It would be of further benefit to vary certain 
parameters while holding others constant, thereby permitting 
conduct of experiments similar to those in the laboratory. 
It is not practicable however, for the traffic engineer to 
study operations under such controlled conditions in the 
field. Therefcre he must resort to methods of analogy, or 
modeling. 

"“Simulaticn prebably provides the most accurate form 

of assessment of a road network, parameters can be 

altered at will, and the resulting changes in delay 

and journey times noted for costing purposes. The 


parameters that can be altered are numerous and 
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include not only traffic signal control policies, but 

also road alterations or the flows and routes taken by 

different categories of vehicles. It is the ease with 
which the network can be altered, controlled and 
assessed that makes simulation so vital to area 

traffic control" (48; p. 464). 

The application of simulaticn techniques enables the study 
of complex traffic systems in the labcratory rather than in 
the field. It is usually faster and less expensive than the 
testing of real system and in many cases enables the _ study 
of system characteristics prior to construction of the 
facility. In the study of traffic simulation, it is 
possible to represent traffic of a particular characteristic 
desired and in the quantities desired, whereas to obtain the 
characteristics in the field may be very difficult. 
Furthermore, there are substantial hazards in conducting 
field trials whereas such situations can be studied by 
Simulation without risk. 

Blackmore [5] at the Road Research Laboratory has 
conducted experiments in the field by setting up 
experimental tracks and briefing the participating drivers 
on how to drive through the tracks. It is evident that the 
success of such a field experiment depends on the co- 
operation of the participating drivers, and that such 
experiments could be very expensive. 


"For the most part it can be said that the goals 
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achievable by simulation in the traffic flow process 
are clearcut and rewarding. Simulation is an ideal 
technique for traffic research. The simulation model 
is not just another means for accomplishing what we 
can do today but it is a tool for solving prcblems 


which cannot be solved today" [12; p. 286]. 


2-4.2 Methods of Simulation in Traffic Flow 

There are two general types of computers namely: 
Analog (continous-variable) and Digital (discrete-variable) 
computers, and both types are employed in simulation 
studies. A simulation technique is classified as analog or 
digital, depending on the way the problem is formulated and 
the type of computer needed for the formulation and 


solution. 


Gat ace it Analog Simulation 

An analog computer is one in which computation is 
performed by varying the state of some physical elements in 
which the variables are continuous. When analog simulation 
is employed, all parts of the system must be simulated 
Simultaneously. The parallel nature of analog simulation 
therefore makes it closer to the traffic behaviour than does 
digital simulation. 

Each component or function of the analog simulation 
system must be simulated by one or more components in the 


computer. This nature of analog simulation always requires 
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the addition of more computer equipment as the system being 
Simulated becomes more and mcre complex. For small systems 
this requirement is not sericus, but for the study of large 
systems the addition of more simulator elements can become 
expensive. Moreover, in many cases further additions of 
Simulator elements may not ke feasible because there is a 
practical limit to the number of simulator elements that 
will work together satisfactorily. 

Furthermore, the accuracy of the analog computer is 
limited to the accuracy of the physical components involved. 
Thus, the use of analog simulation in traffic studies 
becomes attractive when a special part of a general system 
is to be scrutinized, and this leads to the building of a 
special purpose model; but analog technique becomes less 
attractive as the problem under study becomes more and more 


generalized. 


2EUS2S2 Digital Simulation 

Digital simulation is characterized by the use of a 
digital computer. Whereas the analog computer must handle 
all elements of the simulation system simultaneously (in 
parallel) the digital computer handles elements of the 
Simulation system one after another (in series). eto aes 
way, an increase in the system complexity results in an 


increase in the time required for computation. Although 


accuracy may not necessarily be an important factor in the 
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Simulation, it is possible in digital simulation to reach 
any degree of accuracy desired by carrying out more 
computations. 

In analog simulation, the mathematical models used 
must be those which involve differential equations or which 
can be made te look as thcugh they involve differential 
equations. In digital simulation, the analyst has the 
choice as to the representation of traffic on the roadway 
depending on the detail of simulation required. There are 
two major methods of representation in a traffic simulation, 
the microscopic representation and the macroscopic 
representation. For single intersections, the microscopic 
approach of detailed vehicle characteristic assessment is 
usually employed. When considering vast intersection 
networks, the large amount of real traffic and area involved 
in the simulation may obviate the use of the microscopic 
approach. In such situations, the macroscopic approach is 
imperative because of computer time and core storage 
requirements. Again, the choice will be a function cf the 
available facility and the desired level of detail required 


in the model. 
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2.4.3 Vehicle-roadway Representation 

In digital simulaticn, where there is a flow of 
discrete objects such as vehicles on the roadway, there are 
two general ways in which the objects may be represented for 
the purposes cf simulaticn in the computer, the physical 


representation and the memorandum representation. 


2-4.3.1 Physical representation 

With the physical representation, one or more binary 
digits are assigned to represent the presence, position and 
perhaps the size of the item or vehicle to be simulated. 
Certain areas of the computer memory are assigned and 
organized in such a way as to represent the flow network, 
and the groups of binary digits representing the items are 
caused to flow in the network by suitable manipulations. 

The inplementation of the physical method of vehicle 
representation requires considerable familiarity with the 
computer hardware configuration, and as such is not a 
suitable methcd to choose when developing a simulation model 
which is meant to be a tool for non-computer-oriented 


personnel. 
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2.4.3.2 Memorandum representation 

Memorandum representation consists of recording all 
conditions pertaining to a given vehicle where 
representation is by bits, bytes or words of computer memory 
uniquely associated with this vehicle. These defined memory 
locations are usually referred to as vehicle parameters. 
When the vehicle parameters are stored using a single coded 
computer word, the parts of the word signifying different 
attributes of the vehicle<driver characteristics can be 
extracted and interpreted using suitable computer routines. 
This method of attribute extraction and storage to and from 
a computer word (usually referred to as ‘unpacking and 
packing of a word’) could be time consuming, especially when 
certain attributes such as the position and velocity of the 
vehicle are to be updated every scan interval of the 
Simulation. When computer memory preservation is not 
crucial, especially with the advent of virtual memory in 
modern computers, it would be more efficient to let one 
computer word store a single attribute and thereby reduce 
the time used in computation. 

The memorandum method is easier to understand and to 
program for the computer than the physical method. The 
status of each vehicle is kept in the memory circuit of the 
computer when using the memorandum method. The data for 
each vehicle usually include position (distance from 


starting point cr reference point), lane, desired speed, 
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actual speed, turn requirements (points at which turns are 
to be made). Other desirable information would include the 
length or class of vehicle (whether a car or a truck), 
normal acceleration and deceleration, passing 
Characteristics of drivers, and the time at which the 
vehicle entered the system. 

One computer word is usually assigned to each vehicle 
to conserve ccmputer storage space. In the memorandum 
representation, a register cr memory cell of the computer is 
used as a clock counter. The clock counter is advanced at 
each simulation cycle, and the memoranda of vehicles and 
other data are scanned to determine the kinds of actions to 
be taken by different parts of the system. These actions 
can be summarized in the following decision blocks: 

a) Is it time for a new vehicle to enter the system? 
b) Is this vehicle going to encounter another vehicle 
during the next unit of time? 
c) Is it safe to pass the leader? 
d) Is it time for a vehicle to select a routing lane for 
intersection maneuver? 
Once such questions have been answered, each vehicle can be 
advanced by changing the record to show its position and 
speed one time unit later. The position one time unit later 
is calculated by multiplying the vehicle's speed by the 
length of the unit of time and adding the product to the 


present position of vehicle. 


gaiees4 
sd? doidw 28 
sesasiatadionel 
efotiss Hoes of beayiaes cima dee rao =a», 
eupaiidece Sif wT .sdage ‘Spbioz= xesuquo> evassuoD “oF 
ai istuqeeo edt Yo Lies yaowsa xo eretpar’s a0 brstaeas sqa7 
+e Weonsvbe ef aetdvov Avofd ‘sdT <tetauos J20fo » as bean 
boa asloidey to sbusvoRos sit Bas winyo” noitsluwie do89 
cf snobton to 2ictd ott satanetob ot Beans2 one steb ‘zedso 
encifos seer? .meteye oft fo stanq spetsTttb yd deded- od 
:ey00h soiatssh préwolfo?, od? i-hestassaue ed BBO 

Tesiaye edt schiie da eteidos Sant aed outt #2 at (e* 
sisidev aaddens sze¢nvoome oF Pwiop efoideyv abd et od” 
Ssui2 to tinw sxon ont pakseb’ 

trebee! sus. aai6y o+ aise th a 40 

40% snsl poltvox n sopLsa os alpiden 5 t97 ondt ar & 
Crevdenna aoltosetegak =e 

od nan sistdsy dose tszewias doed) oved act zeenp dies wom 
Hite aobstiaoy- Ft wpde 0} bagoss aa poivasd> yd DsoasvBs 
sored sinw abt ano; nakdtegd a “Gta 4 ao ‘oats wae 0K 
oss 46 beege 2 
i ilps att 


40 


The memorandum method may be varied by allocating one 
memory cell to each unit blicck. If a vehicle is in a given 
unit block then data concerning the vehicle may be stored in 
the corresponding memory cell. Thus it may not be necessary 
to otherwise record the position and lane of each individual 
vehicle in the system. In this way, different forms of 
memorandum representation cculd be applied depending on the 
form of distribution of traffic in the network or the 
individual intersection as required by the analyst. 

The form in which the network or the individual 
intersection is represented in a computer program will 
determine the fcllowing: 

1. The amount of information concerning individual 
vehicles that can be stored; for example, position, 
type, speed, entry lane, destination, travelling times 
and so on. 

2. The computing time and storage requirements of the 
program. 

3. The structure of the simulation, taking into account 
the updating routines and routing rules in the systen. 

The above features in computer programs demand some form of 
memorandum representations, typical examples of which are: 

a) Individual vehicle representation, where each vehicle 
is allocated an initial position, speed and 
destinaticn. Its position and speed being upfdated 


each computer cycle. 
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b) Vehicle spaces, where the network is divided into 
vehicle spaces and a number stored in each address 
representing a vehicle space indicates the speed of 
the vehicle occupying that space. Usually a zero 
indicates a vacant space. In this type of 
representaticn, vehicles have no identity and 
intersection routings are based upon Simple 
deterministic or probabilistic rules. 

c) Platoons, where the network is divided into vehicle 
Spaces as in (b); the spaces in a lane are then 
grouped into blocks consisting alternatively of 
vehicles, spaces, vehicles. 

The number of vehicles, or spaces in each block is then 
stored in successive computer addresses. For each lane, the 
occupancy of the first block is stored at another address 
(zero for spaces, unity for vehicles). 

Gerlough and Wagner [19] have applied the concept of 
the individual vehicle representation to the detail 
simulaticn of an individual intersection. Longley [34] has 
also applied the representation of vehicle and space 
platoons to the simulaticn otf signalized intersection 
network aimed at ccntrolling and balancing queues at the 


various approach lanes. 
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2.4.3.2.1 Individual Vehicle Representation: 

Gerlough and Wagner ({19] of the Planning and Research 
Corporation at Los Angeles have developed and implemented a 
detailed microscopic simulation model of an individual 
intersection. Their simulation employed the individual 
vehicle representation method and was implemented with a 
General Purpose Programming Language -- FORTRAN II. The 
Simulation model of Gerlcugh and Wagner is an isclated 
Signalized intersection in which the roadway is represented 
as an orthogonal intersecticn of two six-lane bidirectional 
roadways. 

Each artery constitutes a co-ordinate system whose 

Origin happens to be the intersection entry line, and in 
which the positive scale is the direction of flow from the 
artery in question. All negative co-ordinates are 
considered as part of approach leg and all positive co- 
ordinates are part of an exit leg. 
Upon entry into the system, vehicles have negative 
positions, travel to zero pcsition at the entry line to the 
intersection, and begin a journey on a positive scale as 
they enter the intersection. 

Vehicles are generated by lane at the predefined 
system entry feints according to a composite exponential 
distribution developed by Kell [29]. For a particular 
direction, such as the Northbound direction, each approach 


lane of the artery forms a list with its continuous lane on 
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the exit leg of the Southbcund artery; so that for the six- 
lane bidirectional intersection, the maximum number of lists 
that would be maintained is 12. Thus, each list is made up 
of an approach portion and exit portion with intermediate 
portion being the intersecticn estate. Vehicle processing 
through the system is done by means of pointers to the first 
and last vehicles of approach portions and exit portions of 
the lists. 

The number of vehicles in the approach lanes and the 
exit lanes cf a list at any time are stored in some 
variables. By identifying each vehicle with a lane-list 
corresponding to its directional co-ordinate systen, 
straight through vehicles are simply processed from the 
approach porticn cf a list to the exit portion of the same 
list and turning movements are accomplished by transfering 
vehicles from a particular lane list (originating list) to 
another (destination list). 

Maximum turning velocities for all turn movements are 
derived according to the turn radii of each particular 
turning path, and this velocity is assumed to be sustained 
throughout the turn process when turning can be nade 
unimpeded. Right turns are assumed to follow the same path 
whereas left turns have two configurations--one for free 
flow turns and the other for delayed turns when there is a 
conflict in turning due to cncoming vehicles. Since maximum 


turn velocities are derived based solely on turn radii, 
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there are two types of maximum turns velocities for left- 
turning vehicles: one for free flow left-turners and the 
other for delayed left-turners. 

Vehicle mcticn in the system is evaluated by applying 
some acceleration equations to the individual vehicles. 
These equations are referred to as the ‘'car-following and 
the free-behavicur!' equations. 

A vehicle in the system is either a leader or a 
follower, it is automatically a leader if it is the first 
vehicle on a list, but it could still be considered a leader 
if it is not the first cn a list and still its free flow 
motion is uninhibited by preceding vehicles; otherwise it is 
a follower. During each period of the simulation, a vehicle 
is updated by applying the appropriate equation to it 
according to the position cf the vehicle in relation to the 
rest of the vehicles in the system and to other decision 
models such as signal system. 

Slowing down and or stopping occur in a variety of 
situations such as red light approaches, amber light 
approaches and blocked turns. Upon generation, each vehicle 
is tagged with desired velocity, desired acceleration and 
deceleration; all obtained from precoded distributions whose 
means are supplied as input to the model. 

Queue statistics are provided for in the model and the 
underlying assumption here is that vehicles become members 


of a queue when they stop behind other vehicles which are 
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already in the queue. For a leader of an approach lane of a 
list, it enters into a queue when it is required to stop at 
the intersection stop line due to one of the stcpping 
criteria. A great deal of analytical statistics are 
provided for in the system to be used as parameters for the 
measurement of effectiveness. Included are such things as 
the number of stops, stopped time per vehicle, maximun 
stopped delay, mean stopped delay and so on. 

Mean system delay (referred to by Gerlough and Wagner 
[19] as "primary measures cf effectiveness') was correlated 
to most of the rest of the statistics gathered in the system 
collectively referred to by Gerlough and Wagner as 
"secondary measures of effectiveness’. The seccndary 
measures of effectiveness included the mean stopped delay, 
mean stopped delay per stopped vehicle, mean queue length, 
mean delay in queue, and proportion of vehicles stcpped. 
The results of the regressicn analysis showed that all the 
five secondary measures of effectiveness were very strongly 
correlated with the mean system delay, the primary measure 


of effectiveness. 


2-4.3.2.2 Vehicle and Space Platoons Representation: 
Longley [34] has applied the representation of vehicle 
and space platoons to the Simulation of signalized 
intersection networks aimed at controlling and balancing 
queues in the various approach lanes. In his 


representation, current vehicle speeds are not stored, but 
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are allocated on the basis cf headway only. The apparent 
major objection to this technique as pointed out by Longley 
himself is the fact that vehicle acceleration constraints 
are neglected. But to use it as a readily available tool, 
the traffic engineer does need to make field surveys in 
order to be able to estimate accelerations and decelerations 
of drivers using the roadway, so that the crude model of 
Longley is deemed satisfactcry. 

In the simulation, speeds are restricted to the 
discrete values of 0, 15, or 30 miles per hour only but 
speeds at the intersection are restricted to 15 miles per 
hour only so that acceleration from 0 to 30 miles per hour 
in one second, which would ctherwise be unrealistic, is 
avoided. Longley, by these speed assignments, assumes that 
travel speeds of 0 tc 10 miles per hour are effectively 
stopped; speeds of 11 to 20 miles per hour are equivalent to 
15 miles per hour, and speeds of 21 to 30 miles per hour. 

are equivalent to 30 miles per hour. 

These discrete speeds look unrealistic at first 
glance, but for simulation of large intersection models or 
networks, the above speed assignments approach the 
continuous nature of actual driver speeds based on 
acceleration and deceleraticn rates; yet much computation 
time is saved by by-passing the detail computation of travel 
speeds for individual vehicles based on acceleration and 


deceleration models. However, the speed assignments could 
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approximate closely the continuous speed distribution 
scheme, by widening the range of choice of speeds from 
07 15,30), qbotgsa yg07310 +e 207 S0s0teevenn0y 57 A0, A5S¢e2 04525; 
30 miles per hour and still maintain the time-saving headway 
method of speed assignment (determination) . 

Each vehicle space in lLongley's model represents a 
distance of 22 feet and vehicle positions are updated in 
each computer cycle representing 1 second real time. The 
updating routine is performed in two phases. In the first 
phase all lane vehicles except those queueing at the 
intersection stcp lines, are moved up one vehicle space. In 
the second phase platoons of vehicles with individual 
headways of two or more vehicle spaces are moved along a 
further 22 feet corresponding to speeds of 30 miles per 
hour. With British cars the effective car length of 22 feet 
as chosen by Longley is satisfactory. In North America 
however, a more realistic effective car length would be 25 
feet. 

Movement of vehicles in the four intersections making 
up the network is carried out in two stages: The occupancy 
of the first vehicle space from which the vehicle would 
leave is determined, and a check is made to see if a space 
exists at the rear of the receiving lane in the 
intersection. this first procedure is necessary in order to 
avoid the process of making all necessary checks in 


preparation for a merging vehicle from an approach when in 
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actual fact, no vehicle might have reached the point of 
entering the intersection yet. 

Thus the knowledge of the occupancy of the first space 
is a necessity in Longley's model, just as Gerlough and 
Wagner [19] store the information about the first vehicle in 
any list. 

The intersections are considered to be 66 feet square 
and non-turning traffic mcves across them at a speed of 15 
miles per hour. Vehicles in the intersection square move 
across it when the light is amber or green, and enter the 
center lane of the receiving road, provided that vacant 
Space exists at the rear end. Entry into the intersection 
was restricted by the condition that a vehicle could enter 
the intersection only if there are a total of less than four 
vehicles in the intersection, and four vehicle spaces at the 
rear end of the receiving lane. This rule ensured that no 
vehicle was stranded in the intersection so long as the 
amber time equaled or exceeded 3 seconds; the vehicle could 
always clear the intersecticn during the amber phase at the 
prescribed speed of 15 miles per hour or 22 feet per second. 
This is because Longley does not allow a vehicle to enter 
the intersecticn if the signal goes amber, so that with an 
amber time of at least 3 seconds, the three vehicle spaces 
in the intersection lane would be cleared at the constant 
speed of 15 miles per hour ky any vehicle which happened to 


be caught in the intersecticn before the signal went amber. 
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But the condition that vehicles are never allowed to 
enter the intersections when the signal goes amber is too 
strict and unrealistic, since from the speed distributions 
in his model, vehicles could reach the intersection stop 
line at speed of 30 miles per hour when the signal goes 
amber, and at such a speed it is unrealistic not to allow 
some of the vehicles on the platoon to run the amber as 
happens in real traffic situations. It would seem more 
reasonable to stop vehicles from running the amber when a 
vehicle is already stopped before the signal went amber, but 
that some probabilistic way be employed to let a vehicle run 
an amber if the vehicle is already in motion and at the 


intersection stcp line at the time the signal goes amber. 


2.4.4 Scanning Techniques 

In digital simulaticn, whether by standard general 
purpose programming language or by a general purpose 
simulation language, time is usually divided into discrete 
units, and since it is impossible to examine, within the 
digital computer, all parts of the system simultanecusly, 
some form of scanning must be employed in the simulation. 
Simulaticn packages such as_ GPSS maintain this fcrm of 
clock updating and correlate events automatically. With the 
standard programming languages, the clock updating routine 
is the responsibility of the analyst or the user. 


Whether the system accounts for the clock updating or 
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the user accounts for clock updating, there are two general 
scanning methods available. These are the periodic scan and 
the event scan. 
“An event-scan program is essentially asking ‘what 
happens next?'; whereas the periodic-scan program asks 
‘what will the situation be one time unit from now?!" 


ei2seps (27408¢ 


2-4.4.1 Periodic Scan 

‘Periodic scanning’ consists of scanning and updating 
the entire system cnce during each unit of time. This 
technique is very straightforward and is easy to program. 
For example when periodic scanning is employed in the 
Simulation of an individual intersection, then during each 
particular scan interval every vehicle is updated in turn 
whether or not it is scheduled to perform any action during 
the current scan interval. When all vehicles have _ been 
updated in terms cf accelerating them, decelerating then, 
stopping them or whatever, then all generating points are 
checked for possible entry of new vehicles into the systen. 
In short, the entire simulation system is scanned in any one 
particular interval before the clock is advanced to the next 


Simulation interval. 
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2.4.4.2 Event Scan 

‘Event scanning' consists of determining the next 

event of significance by extrapolation and moving the clock 
to this event without any intervening scans. 
For example, when an event scanning technique is employed in 
simulaticn of traffic circle or any unsignalized 
intersection, vehicles gueveing up behind other vehicles at 
the circle entry line or at the minor roads of the 
unsignalized intersection would not be scanned during each 
scanning interval. Instead, the imminent significant events 
such as the acceptance of gap by a leader of queued vehicles 
and eventual entry into the circle, the entry of new 
vehicles into the system or the exiting of vehicles from the 
circle, would always be determined and stored and the 
earliest one to cccur selected. The occurance of this next 
Significant event may alter the possibility or timing of 
other events that had been listed, so that a new set of 
events and times be calculated. The event-scanning 
technique saves a lot of execution time but requires more 
programming complexities and more computer storage. 

Usually, the periodic-scan and the event-scan methods 
are partially ccmbined in order to produce a program that is 
best suited to the problem at hand. This is true especially 
in the simulation of traffic phenomena at unsignalized 
intersections, uncontrolled intersections, or traffic 


circles. Since vehicles cccupying different parts of the 
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roadway in such a simulation system would have different 
priorities based solely on their positions on the roadway; 
it would be necessary tc update priority-vehicles first 
before the non-priority vehicles. The idea of updating 
priority-vehicles first, especially the leaders in the same 
priority lanes, is ewent-scanning in nature, where the 
earliest imminent events are always selected first for 
updating. On the other hand, the idea of not moving the 
Simulator clock to suit a particular event but instead 
continuing to scan the rest cf the vehicles in the entire 
system after updating the priority-vehicles, is periodic- 
scanning in nature. In simulating traffic phenomena, 
partial combinations of periodic- and event-scanning 
techniques are therefore not only suitable but necessary in 
order to simulate the parallel nature of traffic flow in 


real situations. 


224.5 Other Traffic Simulation Models 

The actual application of computers to the study of 
traffic flow problems was initiated by Goode, Pollmar and 
Wright [20] when in 1956 they developed a model which 
Simulated the movement of traffic through signalized 
intersections. Since the work of Goode et al, traffic 
Simulation studies have expanded considerably with the two 
main methods of approach being the ‘macroscopic method" and 


the 'microscopic method" as already mentioned. 
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2.4.5.1 Macroscopic Simulation Models 

The macrescopic methcd of simulation has been applied 
mainly to large networks of intersections, the sizes of 
which usually obviate the use of microscopic simulation 
methods. The majority of macroscopic traffic simulations 
have the same approach which is the division of the roadway 
into coarse segments and keeping track of the number of 
vehicles in each segment during the entire period of the 
Simulation. Variations of this approach have been used by 
Sakai and Nagca [46], Katz {27] and Gerlough [17]. The 
length of the rcadway segment is dictated by the current 
traffic conditions and the detail desired in the simulation. 
The length of the cell is usually equal to the mean distance 
that an average vehicle would travel in a scan interval, 
which makes the periodic scanning technique suitable for the 
macroscopic simulation of networks. 

Macroscopic network simulation models have varied in 
degree of detail and complexity depending on the primary 
cbjective of the various models. When the modelling of the 
network is used merely asa tool to test the validity or 
effectiveness of some design criteria or proposed strategy, 
then the network is usually simplified to suit the purpose. 
But when the modelling of the network is done solely to 
simulate the kehaviour of an actual traffic network where 
validation of the model would be carried out by comparison 


with real traffic data, as many complexities involved in the 
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real network situation as possible are included in the 
Simulaticn model. 

Examples cf network simulation models in which the 
networks are simplified to study some other design criteria 
are: 

1. The network simulation model by Watjen [51] of the 
Technical University of Denmark. The aim of this study was 
to investigate the delay of vehicles traversing three 
intersections using a variety of synchronizing offsets. The 
secondary objective of his simulation was to study the 
effectiveness of the application of GPSS to traffic 
Simulation studies. Watjen*s model comprised only a single 
lane of one-way traffic and a fixed-cycle signal scheme at 
each of the three linearly arranged intersections. The only 
variables in the model were the offsets of signal two from 
Signal one and that of signal three from signal two. The 
model is of course oversimplified and in fact too simplified 
to be useful; however, the objectives of the study were 
realized. The study shows that GPSS/360 is useful for such 
small scale, non-detailed simulation and that inter-signal 
influences in a network should be considered when simulating 
a large network of signalized intersections. 

2. Similarly, the network employed by Longley [34] in his 
simulation was four symmetrical orthogonal intersections 
which formed a square. In this simulation also, the primary 


objective is the queue balancing strategy, so that the 
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complexity and realism in the network itself is a secondary 
matter. 

However, in other major network simulation studies 
such as that of Katz [27], where model validation is dene by 
comparison with real traffic data from specific locations, 
all complexities and details of the network are included in 
the model, because the primary objective in such models is 
the study of behaviour of vehicles in a real traffic 
network. The advantage in the use of macroscopic 
representations lies in the efficient usage of computer 
storage and the speed of computation, and are applicable 
particularly to large networks where slight errors are 
negligible; however they suffer the disadvantage of lack of 


model realism and detail. 


2-4.5.2 Microscopic Simulation Models 
Microscopic methods are usually employed in the simulation 
of individual intersecticns where detailed realism is 
important. The obvious disadvantage in the use of 
microscopic methods for simulating traffic behaviour lies in 
the inefficient use of computer space and time, but when the 
model is not very big or when computer time and space are 
not crucial, the use of micrcescopic methods should always be 
preferred since microscopic models tend to convey more 
understanding than macroscopic models. The most well- 


defined microscopic simulaticn model is the one by Gerlough 
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and Wagner [19] as described in section 2.4.3.2. 


2.5 Traffic Circles and Traffic Signals 

At high traffic flcws, the most convenient form of 
intersection control is a grade-separated one, since with 
grade-separation only turning traffic encounter any fcrm of 
conflicts, and even this form of conflict is treated by 
imposing a priority intersection such as a 2-way STOP or by 
imposing a traffic circle or a traffic signal. 

Grade separation is employed on freeways and other 
intersections where flows are very heavy, because of the 
high importance attached to the avoidance of any ferm of 
slowing down or delay to main road traffic. The Highway 
Capacity Manual [26] implies that in a rural area at least, 
grade separaticn is justified sige total demand at 
intersection exceeds 4000 vehicles per hour. While grade- 
separated intersections are ideal from the driver's point of 
view, they have major disadvantages such as: 

1. they are extremely expensive to construct; 
2. the amount of land taken up is very high; 
3. interference with pedestrian movement is maximized by 
grade-separated structures. 
The last two disadvantages of grade-separated intersections 
are more pronounced in urban areas with the additional 
disadvantage that grade separation is particularly intrusive 


enviromentally, and may have adverse social effects. 
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In urban areas therefore it is particularly important 
to exploit to the full, the potentialities of at-grade 
intersections. The two najor forms of at-grade 
intersections are Ssignal-controlled intersections and 
traffic circles. The choice between signal-controlled 
intersections and traffic circles is evident from the 
following characteristics of the two forms of control as 
given by Millard [37]: 

a) The land taken up by the ordinary conventional large 
central island traffic circle is greater than that for 
traffic signals. In urban situations, where land 
acquisiticn may be crucial, this is often the deciding 
factor for the choice of signals over circles. 

b) Where flows are unbalanced, particularly if traffic 
entering from one arm greatly exceeds traffic leaving 
by it, as cccur at some intersections during peak 
hours; traffic circles cause undue delays to traffic 
from the minor streets. 

c) Left-turners are as a rule very poorly treated by 
signal-controlled intersections. Usually either a 
special phase is needed or they must queue up in an 
approach. Either procedure delays them and reduces 
the capacity of the intersection. When the number is 
large, a traffic circle treatment is to be preferred. 

d) Three-way intersections are not satisfactorily treated 


by traffic lights, especially if the flows are 
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balanced. The same is true to a lesser degree, at 

five- and more-way intersections. 
The above characteristics spell out the fact that while 
traffic circles may be regarded with disfavour due to 
substantial delays liable tc be experienced by drivers when 
capacity is exceeded, they cannot be ignored or eliminated 
completely from urban street network if utmost efficiency is 
to be achieved in the systen. 

Chapter III describes a model which can be used to 
study the operation and performances of traffic circles of 


different configurations. 


Chapter Til 


Model Developement 


3.1 Initial Consideration 

Most traffic simulation models are claimed to be 
tools, aimed at helping the traffic engineer in the 
evaluation of cne traffic situation or the other. However, 
in the majority of these models, the kind of data required, 
such as the acceleration and deceleration constraints, 
vehicle-driver reaction times and other data of the like, 
make such models impracticable due to difficulty with which 
the data could ke cbtained. It was decided therefore to 
develop a traffic simulaticn model which would meet the kind 
of data readily available to the traffic engineer and be 
accurate enough to give a true indication of the traffic 
Situation. 

To assess the perfcrmance of a traffic circle of 
assumed configuration, the random nature of the distribution 
of turning vehicles within the circle traffic stream should 
be taken into account ard included in the analysis. fhe 
model which was built to meet the above specifications was 
given the acronym TRACISM (IRAffic CIrcle Simulation Model). 
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In TRACISM the traffic engineer need not specify the 
acceleration and deceleraticn of the drivers occupying the 
two-dimensional roadway. At the beginning of each scan 
interval vehicles are deterministically reassigned 
velocities taking into account only their current headways 
and velocities. 

In this model, the general procedure followed to 
obtain a value cf stochastic variable is to transform by an 
appropriate functicn, a random sample from a uniform 
distribution in the interval 0 to 1. These transforming 
functions are in the case of driver desired velocity, a 
normal distribution and for interarrival time gaps, a 
negative exponential distribution (see Appendix A). After 
generating an arrival, the Monte Carlo technique is used to 
assign destinations (turning cr straight through) to the 


vehicles of the arriving traffic stream. 


3.2 Primary Model Criteria 

The objective of this thesis in the area of simulation 
was to develop a model with sufficient accuracy to describe 
the operation of a traffic circle under a number of 
controlled conditions. The following criteria based cn the 
above objective cf the thesis, have been established for the 


model. 
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General Criteria 

The model shculd be flexible enough to allow 
variations of physical, control and operational 
conditions with minimum effort. To this end, the 
number of circle arms, lanes and the general 
geometrics of the circle should be made variables so 
that they can be changed as required. 

The model should realistically and accurately simulate 
System conditions for the purposes of the study being 
undertaken. This will reguire the use of stochastic 
models to describe operational and driver behaviour 
variations at the traffic circle. The stochastic 
models will include the choice of circle routing 
lanes, the choice of exit lanes and the acceptance of 
gaps in the circle traffic stream from the apfroach 
streams. 

The model should exclude refinements which do not add 
significantly to the accuracy of the simulation 
phenomenon. Such refinements as the actual diameter 
of the central island of the circle, and of the 
circulating lanes, which do not add significantly to 
the operational efficiency of the model, should 
therefore be excluded from the model. 

The input should be based on operational patterns and 
Characteristics which can be readily measured and 


checked by the traffic engineer. In this respect, 
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such data as driver reaction times, accelerations and 
decelerations which are not so obviously and readily 
measured by the engineer, should be eliminated if 
possible. 

The output should be designed to allow testing and 
analysis in the simplest possible manner. As a tool 
therefore, the output should contain those statistics 
which are relevant to the engineer; typical of such 
outputs are measures like the time spent in apfroach 
queues, the demand at the circle (number of vehicles 
entering the system), the capacity of the circle 
(number cf vehicles leaving the system within the 
Simulation period), mean system transit times, mean 
system delays and so on. 

The amount of data and its processing for a particular 


investigation should be minimal. 


System Criteria 

The model should be flexible enough to describe 
operation of a traffic circle for any particular 
vehicle, aS an average per vehicle or as a total for 
the system. Thus, the model should be able to 
Simulate trucks only or cars only or the mixture 
thereof. 

The system should be sensitive enough to show 


operational variation with a significant alteration at 
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any one point in the system. For instance, if the 
frequency of arrivals at a particular artery changes 
Significantly during some specific short period of 
time, such a variation in the demand should reflect on 
the delay to vehicles in this particular approach 
direction, the queue lengths in the lanes of this 
approach direction, and in the system as a whole, and 
all the measures of effectiveness which depend on the 


traffic demand at that particular approach direction. 


3.3 Simulation with FORTRAN 
In discrete-event simulation the language used must 
enable the system designer to represent a complex system 
conveniently and comfortably. The representations could 
involve such requirements as solving of equations, 
preserving interrelations, representing decisicn logics and 
physical characteristics cf the system. The languages 
employed vary from Assembler languages through the general 
purpose programming languages (FORTRAN, PL/1, ALGOL) to the 
general Purpose Simulation languages (GPSS, SIMSCRIPT). In 
TRACISM, the general purpose programming language FORTRAN IV 
was chosen as the tool over the available General Purpose 
Simulaticn languages, notably GPSS, for the following 
reasons: 
(a) GPSS/360 is available only on larger machines so that 


the use of such a language would limit the sccpe of 
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the model and hence limit its use as a tool for the 
traffic engineers. 
(b) The speed of executicn in GPSS tends to decrease with 
the growth of model size and complexity; and as a 
tool, it is important that execution be as fast as 
possible to make the use of the model economical. 
(c) It is easier to expand or decrease the size of the 
model through redimensioning of the appropriate 
FORTRAN arrays than it would be for a novice 
programmer to use the REALLOCATE feature in GPSS to 
expand a model programmed in GPSS. 
According to Reitman [45], the general requirements of a 
language in the area of simulation may be reduced to four 
basic characteristics as follows: 

1. Short-term results; 

2. Ability of system to represent real world; 

3. Long-term results; 

4. Effort required. 
FORTRAN meets the above four requirements for the purposes 
of this study as summarized below. 

(a) Short-term results. The system designer must have a 
gcod programming background in the language before he can 
use it for simulation. But FORTRAN, being one of the most 
common languages gives less difficulty as far as background 
in programming is concerned. 


(b) Ability of system to represent real world. Almost 
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any real-world condition could be represented in FORTRAN; 


the effort goes up, however with complexity in a nonlinear 
relationship. 
(i) Logical situations in the model can be represented 
without difficulty. 

(ii) Mathematical capability of the language is excellent. 
There are numerous special-purpose techniques for data 
smoothing, linear programming, and other forms of data 
Manipulations which are both available and accessible 
for the simulation. 

(iii) Maximum model size is completely under the control of 


(Cc) 


the programmer. He can make trade-offs between 
storage hierarchy and speed of execution. 


Long-term results. In this area lies the most 


important advantage of language generality. 


(1) 


(ii) 


(iii) 


Documentation is under the control of the individual. 
There are aids in the form of cross-reference files 
after the individual has thoroughly set up his 
comments. 

Systen designers cther than the original model 
developer can easily follow the logic and detail of 
the simulaticn with less effort. 

Computers of different manufacturers can use the same 
higher-order language program. Usually there is some 
requirement for rework; but in the overall size of 


effort, this would be considered minor. 
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(d) Effort required. Less effort is required to carry 
out the simulation in FORTRAN than it would be for the rest 
of the General-Purpose programming and simulation languages. 
Several programmers can werk on the simulation in parallel 
if the conventions governing the transfer of data between 
subroutines are well planned in advance. In this respect, a 
Simulation program written in FORTRAN, in modular fashion, 
can easily be altered by anyone with a working knowledge of 
FORTRAN, and in particular, parts of the program can be 


altered without affecting the whole program. 


3.4 TRACISM (LRAffic Circle Simulation Model). 

The Traffic Circle Simulation Model has been developed 
as a tool for the traffic engineer in his evaluaticn of 
traffic circle performance in terms of capacity and other 
measures of effectiveness. The model has therefore been 
designed to be as flexible as possible in terms of ease of 
application and some other general terms; and still give 
results accurate enough fcr the purposes of the evaluation 
study being undertaken. A more detailed description of the 
model in terms of programming logic and input data 
organization is given in Appendix C. This present section 
describes the model at the design, geometric and 


representation levels. 
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3.4.1 The Physical System 

The physical system represented by the model is a 
traffic circle of arbitrary central island and up to six 
circle arms (also called ‘arteries') and 1250 feet of 
approach distance. Entry inte the circle is controlled by 
YIELD signs situated at a car length from the circle entry 
lines of the various approaches. The physical dimensions of 
the circle in terms of central island radius, and the radii 
of the circulating lanes are not supplied to the simulator 
as input; instead the following configuration options are 
specified by the user of the model: 

1. The number of arms at the circle. This could te as 
few as three and as many aS six, but in each case, it is 
assumed that the circle distance from one arm to the other 
is the same for all adjacent arms, so that for a 4-arm 
circle for instance, the distance from one arm to the other 
along a particular circulating lane constitutes a quadrant 
of the circle. Each arm of the circle is made up cf an 
approach porticn and an exit portion so that ideally, the 
number of generating points in the system is always equal to 
the number of arms at the circle. However, an arm of the 
Circle could te transformed into a one-way street into or 
out of the circle by merely suppressing the exit cr approach 
portion of that arm during an input conversion phase, 

2. The approach distance of all approach lanes in the 


system is supplied to the simulator as input. This 
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distance, when specified, is the same for all approach lanes 
and sets up the system boundary for the model. The approach 
distance should be a multiple of the effective length of a 
queued vehicle (also specified by the user) so that the 
maximum number cf vehicles that can be queued up at an 
approach lane at one time could easily be calculated at the 
beginning of the simulation period. 
Currently, the distance is set at 1250 feet so that 
with an effective car length of 25 feet there can be a 
maximum cf 50 queued vehicles in any approach lane during 
any particular time of the simulation, and therefore an 
excess of 50 will indicate an intolerable Situation 
(probably a spill-over into an upstream intersection). The 
system can be expanded however, by simply re-dimensioning 
the appropriate arrays in the program; a process that is 
easily accomplished by anyone with a working knowledge of 
Fortran. 
as The circle circulation lanes and the number of lanes 
at the approach and exit pcertions of each arm are set at 2 
each, but again, this figure can be changed to any desired 
value in the input phase. However, it should be mentioned 
that any number of lanes specified in the input phase 
applies to all the three portions of the physical system 
namely: the approach, exit and circle sections. For 
example, if the number of lanes specified is 3, then it 


means that there are 3 circulating lanes around the circle, 
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3 approach lanes and 3 exit lanes at each arm of the circle. 

4. The effective length of a queued vehicle, also 
known as vehicle space in the system, is specified as part 
of physical configuration of the system, since it affects 
the specificaticn of the approach distance of the circle 
arms. 

The various arms and lanes are numbered in an 
anticlockwise direction. For a 2-lane circle configuration 
all inner lanes on all parts of the system are odd-numbered 
lanes, and all outer lanes are even-numbered lanes. The 
system configuraticn depicting the numbering system in the 


model is shown in Figure 3.1. 


3.4.2 Vehicle Representation 

Vehicles are represented in the system by the 
individual vehicle space method as described in this 
section. A vehicle space is equivalent to an effective car 
length as specified in the input phase, so that the various 
sections of the system are considered as being made up of 
discrete vehicle spaces. However, vehicles are not 
physically moved from one vehicle space to the other within 
the same section of the system. Instead, a vehicle in the 
system is always thought of as having occupied a vehicle 
space in a section of the physical system, and movement is 
accomplished by updating the position parameters cf the 


individual vehicles. 
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Figure 3.1 
System configuration showing all section lanes. 
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Vehicles are generated at the various generating 
points cf the circle arms cn an approach basis and assigned 
to a lane of that approach depending on the predetermined 
destinations. Upon generation, each vehicle has 11 
parameters associated with it. Some of these parameters are 
given initial values upon entry into the system while cthers 
are given values during the course of the vehicle's passage 
through the system and at the time of exit from the systen. 
The parameters that are initialized on entry into the system 
are: the originating approach lane, destination artery, time 
of entry into the system, vehicle desired velocity and the 
current velocity and position of the vehicle. 

During the passage of the vehicle through the systen, 
the current velocity and position are updated continuously 
while cther parameters such as the time of entry into queue 
and the time of entry into the circle are recorded and 
stored in the appropriate parameters of the vehicle. The 
artery by which the vehicle enters the system is not stored 
as a parameter attached to the vehicle, since a knowledge of 
the approach lane number, automatically gives an indication 
of the originating artery of the circle. A vehicle covers 
the approach distance at a speed of 30 miles per hour (the 
maximum allowable speed on the approach portion cf the 
system) unless it follows a slower vehicle or unless it has 
to join a queue at the appreach. Upon entry into the circle 


section of the system, vehicles travel at 20 miles per hour 
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(the maximum speed on the circle) unless they follow slower 
vehicles. When a vehicle comes to the desired exiting 
artery of the circle, it exits from the system by travelling 
one vehicle space along the exit lane. After travelling 
that distance along the exit lane, the time of departure 
from the system is recorded, the rest of the desired 
statistics about the vehicle are gathered by accessing the 
parameters associated with the vehicle, and the vehicle is 
eventually removed from the system. The statistics gathered 
about the individual vehicles are analyzed on per lane, per 
approach and fer artery basis to produce the desired 


measures of effectiveness as described in section 3.4.7. 


3.4.3 Distance Headway Measurement 

The determination of distance headways of vehicles in 
the system is the same for all vehicles on the approach, 
circle and exit portions of the system. When the distance 
headway has been determined the assignment of subsequent 
speeds is different for vehicles in the approach lanes and 
in the circle lanes, because the model assumes different 
speed limits for vehicles on the approach lanes and the 
circle lanes. Presently, the speed limits have been set at 
30 miles per hour on the approach and 20 miles per hour on 
the circle. In distance headway measurements we scan all 
lists in the system, computing Spacing between all vehicles 


during each scan interval. 
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Suppose the average length of a vehicle is 17 feet, 
and the effective length of a queued vehicle is taken to be 
25 feet in the system; this implies that the length of a 
vehicle space is 25 feet and therefore a minimum spacing of 
8 feet between successive vehicles is inherent in the 
system. There are three ways in which the vehicle position 
is reckoned in simulation systems. The position of a 
vehicle in a simulation system could be taken as the 
position concurrent with the front bumper of the vehicle, or 
the rear bumper or even the midsection of the vehicle. With 
any of the above representations a minimum safe distance of 
8 feet between successive vehicles becomes inherent. 

In this model, the positicn of the vehicle at any time 
in the simulation is taken to be the position concurrent 
with the Front Eumper of the vehicle. 

The following procedure of determining distance 
headways of vehicles in the system was adopted: 

Suppose the position of vehicle I, the leader, at time t is 
X feet, that is 

POS(1,t) = X 
and POS(I+1,t) = X* = Position of vehicle I+1, the follower 
at time t. Then the headway H between vehicle I and vehicle 


I+1 during the scan interval at time t is given by 


H = POS(I,t) -(POS(I+1,t) +ELENTA) 


X-(X'+ELENTA) where 


ELENTA is the specified effective length of a queued 
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vehicle, also taken to be the length of a vehicle space in 
the system. The above procedure of calculating headways of 
vehicles in the system always maintains the existence of a 
safe minimum inter-vehicular spacing in the model (as 


illustrated in Figure 3.2). 


3.4.4 Vehicle Speed Determination 

Current vehicle speed and position are continuously 
stored during the simulation period, so that new speeds and 
positions can be computed for the vehicles during the next 
scan interval. 

Speed distributions in the system are approximated as 


summarized below: 


SPEED CODE SPO SE1 SP2 SP3 
VALUE 0 1 2 3 
M.P.H 0 10 20 30 
FEET/SEC 0 15 30 45 


During a scan interval (equivalent to 1 second real 
time) vehicles in the system are either stopped OL 
travelling at 15, 30, or 45 feet per second depending cn the 
available distance headway and the section of the physical 
system being occupied by the vehicles, Around the circle 
however, vehicles are either stopped, or travelling at 15 
feet per second or 30 feet per second (corresponding to 
speeds of 0, 10, or 20 miles per hour respectively); since a 
speed of 30 miles per hour is not permitted around the 


circle. 
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Figure 3.2 
Vehicle positions. 
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As soon as a vehicle stops in the approach lane (that 
is7 when the vehicle attains a speed of SP0), its 
progression is monitored until it is able to merge with the 
circle traffic stream, This accounts for the delay on the 
approach. Suppose that the distance headway H of a vehicle 
has been determined (as in section 3.4.3), then the speeds 
and positions are updated as follows: 

For the leader of an approach lane the updating logic is: 

(a) H < 15 feet VEL(1,t+1) <---SPO 
FOS(1,t+1) <---LENT 

where LENT is the approach distance to be covered. That is, 
if the leader of the approach lane is within 15 feet of the 
circle entry line and there is no acceptable gap in the 
circle traffic stream, then the vehicle is decelerated to a 
stop at the circle stop line during the next time interval 
and awaits the next acceptable gap. However, within the 
distance of 15 feet of the circle, if there is an acceptable 
lag or gap in the circle traffic stream the vehicle would 
accelerate into the circle. 

(b) 15 < H $ 30 feet VEL (1W7tt 17 "<-==SP4q 

POS(1,t+1) <---POS(1,t+1) + 15 
That is, if the distance headway is greater than 15 feet but 
less than 30 feet (if there is a distance headway of at 
least one car length) the vehicle would accelerate or 
decelerate in order to cover a distance of at most 15 feet 


during the next time interval, depending on the current 
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speed of the vehicle. Such an updating process may seem 
unreasonable, but remembering the fact that the model is 
structured towards peak hour flows, it is reasonable to 
assume that queued vehicles move up slowly at a circle 
whenever there is a headway of at least half a car length. 

(c) 30 < H < 75 feet VEL(1,t+1) <---SP2 

POS (1,t+1) <---POS(1,t+1) + 30 
For a distance headway of at least a car length and at most 
three car lengths, a vehicle travels at a speed of 20 miles 
per hour and covers a distance of 30 feet in the next time 
interval. 

(4d) H> 75 feet VEL (1,t#1) <---SP3 

POS(1,t+1) <---POS(1,t+1) + 45 
For a distance headway of more than three car lengths, there 
is enough room for a vehicle tec accelerate to a speed of 30 
Miles per hour, providing that it is not accelerating from 
rest (that is if its current speed is not zero). 

For subsequent vehicles, the determination of vehicle 
speeds and eventual updating of vehicle positions is the 
same as that for the leader of the lane except that no 
vehicle other than the leader is allowed to enter the circle 
stream within the same scan interval of 1 real _ second. 
Vehicle speed and position updating are the same fcr the 
Circle lanes as for the apfroach lanes except that the 
maximum speed around the circle is set at 30 feet per second 


(20 miles per hour) as ofposed to the maximum speed of 45 
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feet per second (30 miles per hour) on the approach. 


3.4.5 Model Processing Order 

At a signal-control intersection there is less 
interplay between drivers, especially when the intersection 
is under a 4-way phasing scheme where every intersection 
routing is prctected. Even in the ordinary 2-way phasing 
scheme, there is less interplay between drivers than at a 
traffic circle. This is because, during the green phase, 
straight through and fright-turning vehicles from the 
Cpposing directions have scme kind of synchronized actions 
of intersection routing without much interference other than 
that due to pedestrian; only left-turning vehicles have any 
form of interplay with opposing oncoming vehicles. Thus, in 
simulating a signal-controlled intersection, all lanes of a 
particular appreach could be processed or updated tefore 
considering ancther set of lanes, and still be within 
reasonable accuracy with this fcrm of serial processing. 

However, at a traffic circle where there is almost no 
evidence of synchronized driver reactions from the 
intersecting traffic streams, parallel processing of 
vehicles at the circle is imperative if any form of accuracy 
is expected in terms of approximation to real situations. 
Since the presence of other vehicles in the circle traffic 
stream, the arrival of other vehicles in the circle and the 


exiting of yet cther vehicles from the circle do dictate 
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what other drivers in the approach lanes do, it is important 
to co-ordinate all mcvements in the systen. 

In this model therefore, the sequencing of simulated 
vehicle maneuvers had to ke performed with some care, since 
an essentially parallel operation in real traffic system is 
Simulated by a serial computation in the digital computer. 
As a general rule, vehicles are only allowed to enter a lane 
after the vehicles in that lane have been updated. 

In order to implement this condition, the exit lanes were 
the first to be updated, and the total sequence employed is 
illustrated in Figure 3.3, and summarized below: 

(i) Each individual maneuver (such as the ‘Inner Approach 
Lanes!) was updated in all four arteries before the 
next set cf lanes was processed. 

(ii) Traffic was considered to have left the system as soon 
as a vehicle had passed one vehicle space along the 
exit lane, so that no computer time was wasted in 
updating outward-bound traffic on the exit roads. 

(iii) The circle arms and rcadways were numbered so as to 
retain uniformity of circle operation in the system. 

(iv) The simulation included the assumption that the exit 
roads were clear so that traffic was always free to 
leave the system. This assumption was guaranteed 
valid by the fact that exit lanes were the first to be 


processed at the beginning of each scan interval. 
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Processing order. 
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3.4.6 Circle Entry 

Approach vehicles at the circle entry line which are 
ready to merge with the circle traffic stream have to 
evaluate the gars in the circle stream by means of some 
deterministic rules. The circulating lanes in the circle 
are divided intc sections, forming separate lanes in terms 
of circle configuration and processing. For a 4-arm circle 
with 2 circulating lanes, division of lanes will result in 8 
separate circle lane sections; for a 5-arm circle with 2 
circulating lanes division of the lanes will result in 10 
separate circle lane sections, and so on. 

Vehicles ready to enter the circle find themselves 
confronted with two streams of traffic to their left and 
right. A potential merger attempts an entry into the circle 
after completing the following safety checks: 

1. that the pertion of the circle where he intends to 

drive tc, is free of vehicles; 

2. that he could drive to his destination in the 

circle before a vehicle from a left lane of the circle 

occupies the space. 
In short, a potential merger from an approach lane always 
defines some form of a safety zone which’ should be 
completely free of vehicles before an entry into the circle 
would be attempted. Different potential mergers define 
different unique safety zones based on their own judgements 


taking into account the abilities of their vehicles and 
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their own responses. This type of decision is usually 
carried out by means of a stochastic model which randomly 
assigns safety zones to all vehicles which are ready to 
merge with the circle traffic stream. But in this model, a 
uniform deterministic way of safety zone evaluation is 
adopted after field observaticns of gap acceptance maneuvers 
at some traffic circles. Measurements were made at some 
specific traffic circles in the City of Edmonton as 
described in Appendix E, and upon those measurements a 
safety zone was selected for all vehicles ready to merge 


with the circle traffic stream. 


3.4.7 Gap and Lag Acceptance 

The acceptance or rejection of a time or space gap is 
a binomial response and is dependent on the size of the gap. 
The minimum time or space gap that a driver accepts is fixed 
for that driver. He will reject all gaps smaller than that 
time interval and accept all gaps larger than that time or 
space gap. There is an evidence that this minimun 
acceptance time or space gap would decrease with time 
pressure and the number cf vehicles in the circle traffic 
strean. 

A lag at the circle may be defined as the time 
interval between the arrival of an approach street vehicle, 
and the arrival thereafter of the first vehicle in the 


circle section at a reference point. A gap at the circle is 
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defined as each time or space headway formed by successive 
crossing of a reference line by circle lane vehicles. 

If the approach street vehicle moves into the circle 
before the arrival of the first circle lane vehicle, the 
driver of the approach street vehicle is said to "accept" 
the lag. If he remains until after the first vehicle 
passes, he has "rejected" the lag. After rejecting a lag, 
he then evaluates the gaps between the successive circle 
lane vehicles. Each gap that he fails to accomodate his 
vehicle into is said to be rejected. The gap that the 
driver finally moves his vehicle into is said to be 


accepted. 


3.4.8 Measures cf Effectiveness 

Measures of effectiveness, sometimes referred to as 
figures of merit, are used in most traffic simulation models 
for both validation and research purposes. For validation 
studies, measures that are readily obtainable in supporting 
field studies are chosen for determination of equivalence 
between the simulated system and actual system performance. 
On the other hand, for research purposes, measures are 
chosen that indicate efficiency and uniformity of traffic 
flow. Such measures usually include spot speeds, average 
acceleration and deceleration, system entry speeds and other 
measures which are cften irrelevant as far as comparison 


with actual data is ccncerned. 
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In designing this model, several different measures 
were considered, most of which are for validation studies 
more than for research purpeses. The measures include the 
following: 

1. directional flows at the circle; 
2. average waiting time in a queue; 
3. mean transit time by route; 

4. mean delay by route; 

5. mean system delay; 

6. mean delay by apfroach; 

7. j%mmean delay by lane; 

8. average queue length 

9. maximum queue lengths; 

10. mean system transit time. 

The measures cf effectiveness indicated above are figures 
that are readily obtained by field counts and surveys and 
some data reduction and analysis, so that validation cf the 
model realism by compariscn of simulated phenomenon with 
real traffic data is easily accomplished. Various tables 
are provided for on lane, approach and artery basis, so that 
any form of queue balancing or equalization of delay 
policies can easily be implemented. Some of the tables are 
also provided for on origin-destination (route) basis so as 
to give an indication cf the effect of circle routing on the 
delays and travel times. 


Using this traffic circle model, a number of 
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simulation runs were carried out to indicate the degree of 
detail required in the sinulation of traffic circles. The 
results of these runs for some typical traffic circle 
configurations are given in the next Chapter together with 


the appropriate refinements to the model. 


3.4.9 Theory cf Control Parameters 
Let (a) be the vector which describes the 
probabilities of vehicle arrival at the various apprcaches 


of the circle arteries. Then (a) is given by 


(a)e= : 


Where n is the number of arteries at the traffic circle. 


Let [t] be the matrix which describes the exiting 
probabilities (also called the turning probabilities) cf the 


arriving vehicles; then the matrix is represented as: 


0 t> se e e ce 
{t) = : ‘ 
tn tha - - 0 


[t] is an n n square matrix with zero diagonal elements, 
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where n is as defined above. 


£ 


7 probability that a vehicle from the approach leg of 
artery i will exit from the circle by the exit leg 
of the j-th artery 


i,j = Vile e eeseeceaendecees n a1 # 5 


i) 
© 
}- 
rh 


tj i= 5 e 


Basing a measure of performance of the traffic circle cn the 
queue lengths and the waiting times in the queues as well as 
the delays experienced by the vehicles during their passage 
through the system; the following queue, waiting time and 
delay matrices are maintained in the model: 

(Q] is a matrix which ccntains the average queue lengths 
maintained at the approach janes. FOG Ss cueecitclie Of 
approach lanes at each of the approach legs, [Q] is an nx 2 


matrix for an n-artery traffic circle, given by 


Q, po 

Q2; Qo2 
[ 0 J = e e 

Qn Qn2 


Similarly, [W] the matrix of waiting times in the approach 


queues is represented as follows: 
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Wi Wo 
Wo, Waa 
Cw) = : i 
Way Boa 
Q; = average length of the j-th queue at the approach 
leqeotearteryoi. 
W; = average waiting time in the j-th queue cf 


the i-th artery. 
(D] is a matrix which describes the vehicle transit 


delays in the system, and is represented by: 


0 Dio se e e D,, 
CD)s= : : 
D., D> e es 0 
Dj .= the average delay (slowing and stopping delays) 


to the vehicle which enters the system by the approach leg 
of artery i and exits from the system by the exit leg of 


artery j.- 


The following relationships based on the above ccntrol 


parameters have been established in the model: 
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2. [W] 
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the average queue lengths and the average waiting 


times in the queues are both functions of the arrival 
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probabilities and exiting probabilities (together 

referred to as ‘generating characteristics’). 

3+ _ (Dd subt((a)7thtde) 

the delay to the vehicles is also a function of the 

generating characteristics. 

Thus varying (a) and [t] would result in varying values of 
{0}, (W] and [D]. It is possible therefore to determine the 
values of (a) and [t] by manipulating the total and turning 
volumes of traffic at the circle that would yield saturation 
values for [Q], [W] and {Dj}, assuming that some limits have 
been placed on [0], [W] and {D] beyond which the use of 
traffic circle would be regarded with disfavour. This kind 
Ofeecircie evaluation study is necessary for the 
implementation of some conversion policies from the use of 
an existing traffic circle to traffic signals or the 
decision to ccnstruct new traffic circles within some 
localities. 

When flows are very light it suffices to control the 
intersection pcints with YIELD Signs, 2-way STOP signs or 
even Y-way STOP signs, depending on sone other 
characteristics of the intersection such as the range of 
visibility from the various approach streams. As flows at 
the intersection get heavy, it becomes necessary to gc into 
more restrictive ccntrol such as the use of traffic circle 
OletractscC 9s1 gals. Lt 1S 2Hporcanteamtnateernreetratcic 


engineer establishes the unsuitability of the use of a 
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traffic circle at the particular intersection before going 
into the installation of traffic signals. 

The suitability or unsuitability of a traffic circle 
with respect to the above performance measures can be 
established by controlling the values of (a) and [t]} to 
study the distribution of [C]) , [W] and [D] . By the use of 
a model with signal control, it will be possible to cbtain 
the distributions of [0] , [W] and [D] under signal control, 
and thereafter establish criteria for the implementaticn of 
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Model Validation 


Validation in traffic flow simulation studies is the 
process whereby the simulaticn model is evaluated by either 
a straight comparison with actual traffic behaviour or 
through some statistical methods, to determine whether the 
model satisfactorily duplicates real traffic behaviour. 
Since it is not the goal in traffic simulation to reproduce 
all minute details in real system, it is necessary to 
establish in the beginning of the Simulation those 
characteristics of real traffic which the model must 
duplicate in order to be considered as a useful model; in 
other werds, which criteria are to be used in validation 
process. 

Theoretically, the model should duplicate the 
characteristics that the traffic engineer uses as design 
criteria or the characteristics that the engineer uses as 
operational criteria. However, the choice of the 
characteristics for validation studies is more often than 
not dictated by the feasibility of their measurements in the 
field, since any traffic flow simulation validation is 
meaningless unless there is a measure of comparison from 
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actual traffic situations. 

The lack of adequate data on the operational measures 
of a traffic circle limited the validation studies tc the 
flow characteristics at the circle and the travel times 
during peak flows. From the analysis of the observed data 
as obtained from the Traffic Division of the City of 
Edmonton Engineering and MTransportation Department, the 
directional flows at the study circles were chosen to he the 
characteristics for operational criteria and hence for 
validation precess. Additional field studies were 
undertaken at specific traffic circles to determine the real 
travel times for some particular origin-destination routes. 
The validation process involved two sequential steps: the 
refinement of the model during test runs; and the comparison 


of performance between the simulated and real system. 


4,1 Testing and Refinement 
The testing of the model was done in two stages: 

First, a test run was made on the simulation model for a 
given set of conditions in the system. The realism of the 
model was then tested in a general way by observing whether 
cr not the outputs were reasonable. Second, a more exacting 
test was applied by simulating conditions for actual 
locations (see Appendix E for the specific locations). The 
Simulation was carried out fcr four different cases of peak- 


hour traffic input conditions as follows: 
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Case 

om Constant Incoming flcw rate, and constant turning 
rates 

C2 Varying Incoming flow rate, and constant turning rates 

C3 Constant Incoming flow rate, and varying turning rates 

c4 Varying Incoming flow rate, and varying turning rates 


For each peak-hour traffic data obtained for a particular 
location, the model was emplcyed for successive runs under 
the conditions as specified by the four cases above. The 
output of a typical run for a specific input data is shown 


in figure 4.1 


The average travel times and the number of vehicles through, 
are catagorized by origin-destination, with the approach 
lane as the origin and the exit artery as the destination. 
In this way the 16 allowable routes in the system are given 
by the 16 origin-destination codes starting from 1 through 
16. The allowable routes are the right-turn maneuvers by 
cuter lane vehicles, the left-turn maneuvers by inner lane 
vehicles only and straight-through maneuvers by both inner 


and outer lane vehicles. 
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Run number 


Duration of run 


3600 secs. 


Simulation time interval 1 secs. 


Current time of simulation 2700 secs. 


Fill time 


900 secs. 


Average Travel Times/Vehicles through 


OD-Code 
1 
71 
3 


15 


16 


Time (secs) 
20 
26 
19 
26 
25 
19 
15 
21 
26 
23 
84 
43 
19 
25 
19 


23 


Figure 4.1 


Sample Program Output 
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Figure 4.1 (ccntd.) 


Avg. Queue Length/Max. Queue Length/Avg. Wait in Queue 


Lane Avg. Que. Lnth. Max. Que. Lnth. Avg. Wait 
1 1 6 8 
2 1 7 40 
3 1 5 21 
u 1 6 19 
5 2 13 30 
6 1 15 75 
7 1 5 16 
8 1 9 21 


Vehicle entry/Turn maneuvers/Output flows 


Artery Entry L£ft. Thru Rgt. Exit Flow 


1 hs 354 347 73 763 
2 715 110 581 28 956 
3 B4E 266 S44 37 595 
4 744 20 617 109 Sa2 


Mean System travel time = 26 seconds 


Mean System delay = 6 seconds 


Figure 4.1 


Sample Program Output 
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For example: 

OD-1 South-north maneuver via Inner lanes 

QD=2 South-West maneuver 

OD-3 South-East maneuver 

OD-4 South-North maneuver via outer lane and so on. 
The queue statistics are gathered on per approach lane 
basis, sc that the outpfut shows the queue statistics for the 
8 approach lanes of a 4-legged 2-lane traffic circle. 

The vehicle entry and exit flows are both recorded on 
an artery basis to indicate how many vehicles enter the 
system by an artery and how many leave by that same artery. 
However, the vehicle entry statistics are further stratified 
into left-, thrcugh-, and right-turn maneuvers. Obviously 
the results are specific to the particular traffic circle 
configuration. However, on the basis of these restricted 
studies the model has demonstrated its ability and 
usefulness as a tool to the traffic engineer. 

The accuracy of the model was tested by comparing the 
outputs of the simulation with field measurements of the 
same parameters. This led to an iterative process whereby 
the model was made successively more accurate through a 
series of refinements. The refinements of the model 
included both parameter tuning and logic modifications. 

For parameter tuning, the effective vehicle length was 
varied between 20 and 30 feet while holding the rest cf the 


system parameters constant. Examination of the various 
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output results and comparison with observed data indicated 
that an effective vehicle length of 25 feet was most 
appropriate. In the logic modification stage, the allowable 
routes in the system were initially made to include U-turn 
Maneuvers at the circle; but the inconsistency of the 
Simulation results with observed data (which was collected 
and analyzed on the assumption of non-U-turners) 
neccesitated the exclusicn of U-turn maneuvering from the 
logic of the system. 

The refinement phase was repeated until results were 
deemed satisfactory. Upon completion of the refinement 
phase the model was employed to investigate the comparison 
of performance between the simulated and real systems. 
Comparison of the results was categorized by origin- 
destination on both the simulated and real systems. 

Table 4.1 shows the comparisons employed for the 
various cases of peak-hour traffic demands. The Table lists 
the computed inflow of traffic and the computed results of 
circle capacity in vehicles per hour. Real data were 
derived from traffic counts and surveys conducted by the 
City of Edmonton Engineering and Transportation Department, 
Traffic Division. The Capacity/Demand Ratio is simply the 
ratio of computed capacity and computed incoming flow, both 
in units of vehicles fer hour (VPH). The reascnable 
agreement shown in the computed and observed values supports 


the validity of the model. 
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Takle 4.1 


Comparison of Simulated and Observed phenomena 


Case cA C2 C3 c4 


Morning Peak 


Computed Incoming Flow, VPH 4097 3807 4106 4121 
Observed Incoming Flow, VPH 3456 3456 3456 3456 
Circle Capacity } Computed 4102 3823 4207 4112 

(VPH) } Observed 3456 3456 3456 3456 
Computed Capacity/Demand Ratio 1.001 1.003 1.020 0.997 


The computed values shown in Table 4.1 were obtained for 4 
different cases of peak hour conditions because the specific 
case (under which the available data was collected) was not 
known. It was not possible to duplicate in the simulation, 
the exact traffic conditions of peak-hour flow during which 
the data was collected; sc periodic flow and turn rates 
supplied to the simulator as inputs were based on estimates 
taken from 15-minute counts of peak-hour traffic flow. 

It is interesting to note that more vehicles left the 
system within the one hour peak period than entered it, as 
indicated by some of the Capacity/Demand Ratios. The 
explanation of this discrepancy is the fact that, for any 
one hour period simulated, a 15-minute period of fill-up was 


simulated prior to that one-hour simulation period. 
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Therefore, when the simulation is started, there are already 
some vehicles in the system. Thus the number that would 
leave the system may or may not reflect the number that 
entered it within the one-hour duration. This fill-up 
concept is very resonable, since it would be unrealistic to 
assume that there would be no vehicles in the system before 
the beginning of a one-hour peak period. Table 4.2 lists 
some of the performance measures computed in the model for 
the four different cases of variable generating 


characteristics for a typical traffic circle configuration. 


The model was not validated for all the measures of 
performance due to the lack of adequate data, but since 
driver attitude indicate that delay is a prime concern, the 
statistics from the system are worth considering for 
efficient use of the traffic circle simulation model for 


more detailed evaluation studies on the traffic circle. 


4.2 Statistical Validation 

The validity of the mcdel has further been examined by 
applying a suitable two-sample statistical testing procedure 
to the system statistics for the data set pairs (simulated 
and observed) obtained. The data set pair for the 
statistical validation studies was the simulated and 


observed average travel times by route. 
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Takle 4.2 


Performance measures 


Case ci ¢c2 ¢3 C4 
Avge System Travel Time (secs) 24 23 24 pM | 
Avg. System Delay (secs) 6 10 5 11 
Max. Queue Length, Northbound (veh) 7 9 6 5 
Max. Queue Length, Westbound (veh) 5 7 6 3) 
Max. Queue Length, Southbound (veh) 14 25 10 8 
Max. Queue Length, Eastbcund (veh) 7 10 9 19 


Avg. Wait in Queue, Northbound (secs) 20 18 14 1S: 
Avg. Wait in Queue, Westbound (secs) 15 17 11 11 
Avg. Wait in queue, Southbcund (secs) 44 62 Ze 16 


Avg. Wait in Queue, Eastbound (secs) 19 28 27 40 


data was collected on one of the study circles for the 
afternoon peak-hour flow, by merely timing the individual 
vehicles from the instant they enter the system (after 
crossing a predetermined system boundary) until they exit 
from the system. The timing was done by an ordinary stop 
watch, so that the timing process could be repeated over and 
over again to make the data collection process attain some 
degree of reasonable accuracy. 

Lack of vantage point for the whole of the approach 


legs and exit legs of the circle limited the field studies 
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to vehicles from one approach leg only. This accounted for 
4 origin-destinaticn routes for the right-turning maneuver, 
left-turning maneuver and straight-through maneuvers from 
either of the approach lanes. However, greater faith was 
placed in the average travel times of the right-turning 
Maneuvers, since visibility was better only for adjacent 


arteries of the circle. 


Let mm; be the ere travel time for the i-th run, i= 
Waa Shy aes 8, observed in the real world fora 
particular origin-destination pair in the model. 

Let u; be the corresponding true, but unknown average 
travel time for the simulation. 

Let ti§ be the average travel time measured for a complete 


set of simulaticn runs. 


The null hypothesis: 
f= Fue, SO 123, fe 8 
is then tested by considering the statistic 
n; = a ew 1 A 7S ple dials on that # is; by 
comparing the real average travel times and the 


Simulated average travel times. 


Now, consider the statistics m- t' = Dj . Application of 
the Central Linit Thecrem implies symmetry of the 
distribution of the statistics Tf eee ene Also, these 
statistics can be assumed to be independent. Again, the 


numbers D; ee ema che may be assumed to be continuous 
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random variables. Hence the Wilcoxon Signed Rank Test may 
be used on the numbers 
De eit, SEL \eidp she ctea 6 
to test the following hypothesis: 
Hoes Mea Us 1—l,2,5¢ tees es cg o 
fom — | 72757 (ose eteO 
The steps involved in the Wilcoxon Signed Rank Test are 


summarized in Table 4.3 


Table 4.3 
Wilcoxcn Signed Rank Test for 
Matched System Travel Time Pairs for OD-3 


RUN NUMBER 
Catches ua Caine ee aati ane ek a ae Cae as Lee. mame mak ee 
{ etl A hh) coe eed on ord Mawes dl Sif meal ot Poe 


a 
{Real Avg. Travel | { | { { { { { | 
{Time (secs) (eS (eS Wii See Ceo ieee ie OG Of 
ae ee Oe bn ae ae es ees eee ee 
{Simulated Avg. { | | | | { { | | 
{Travel Time (secs) { 29] 24{ 31]{ 24{ 29{ 304 27{ 30 
a ee een og pe nny ferme a ee fe Fee Pe ee, 


{Difference Nese oe eee Oe (eee Clee) eee 
p-----------—- se oe ee 
(SlomeoLr pitterence| + |=". {>> {+ "| asta Neneh { 
f+ --4- --+- + --+- + —+-— + -H1 
(Rankeotabitrerence| 691-5" |12° 1-3 4 pests Set Pit 9 | 
f--- $f tt ttt 
{Signed Rank HS fe ge lA Eo LP trek Pca Ua { 
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From Table 4.3 for an effective sample size of 6, the two- 
tailed significance probability corresponding to y= 7 (the 
Smaller sum) is given by 
P = (2) (0.281) = 0.562 
Lor OD-—3e(Ssouth-Frastetrartic) . 
We could also define the test statistic T as the sum of the 


assigned ranks of the positive differences, that is, 


8 
T = DR 
ral 


where R; = 0 if D) is negative 


Rj = the rank assigned if D; is positive 
so that T= )+ =14 


NOW, LOL ec twO-ta1.cd test and a “critical” region) of (size 
<= .05, we accept the null hypothesis H if Tis between 
Y,,, (=1, obtained from Table of Quantiles for Wilcoxcn Signed 


Rank Test) and w, (=20, ottained from Table of Quantiles 


7 
for Wilcoxon Signed Rank Test), or equal to either quantile. 
And since T=14 for the test in Table 4.3 the null hypothesis 
is readily accepted. Hence the hypothesis that the 
Simulated average travel times are equal to the real travel 
times cannot be rejected for OD-3 (the South-East traffic). 
For OD-1 (South-West traffic via Inner lane), OD-2 (South- 
West traffic) and OD-4 (South-North traffic via the Outer 


lane), results were less emphatic, because of lack of 


accurate travel time measurements in the field, but at the 
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S-percent significance level, the hypothesis could still not 
ke rejected. 

On the basis of these four tests, corresponding to the 
various origin-destination routes from an approach it is 
concluded that the model is validated satisfactorily, since 


the logic is the same for all approach legs to the circle. 
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Chapter V 


Conclusions and Recommendations 


5.1 Conclusions on Operational Study 

The feasibility of developing a digital ccmputer 
Simulaticn mecdel of a traffic circle for use as a flexible 
and reliable tcol for determining traffic circle 
performance, has been demonstrated. Field validation of 
some test configurations indicates that faithful results can 
be expected from the use of the model. 

The Simulaticn model has the capacity of accommodating 
variations in the assigned values of all significant design 
and operational parameters of a traffic circle. This 
includes geometric parameters, such as number of lanes, 
approach distances, lengths cf circle lane sections, number 
of circle arms, and traffic demands including vehicle-driver 
characteristics. With this model, rational and reliable 
Operation decisions will be possible under controlled 
conditions, heretofore unoktainable when performing such 


studies at real traffic circle locations. 


5.2 Conclusions on Study Methods 
The results obtained from the use of the model on test 
configurations indicate that although the methods used for 
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the field studies are not all that sophisticated and perhaps 
unacceptable in terms cf accuracy and precision, the 
iterative procedure of simulation refinement studies enabled 
reasonable results to be obtained. The methods used for the 
gap measurements and headway measurements cost almost 
nothing and can therefore be repeated over and over again as 
required in the simulation. However, to provide a permanent 
study record and to facilitate desired exactness of 
measurements, time-lapse photography should be used as the 
most appropriate means of recording gaps, lags and headways. 

The flexibility of the gap/lag acceptance routine has 
been built into the model so that the gap acceptance logic 
can be changed to suit more faithful field studies that may 
be undertaken in future. The present logic in the systen 
employs the simple idea of a vehicle accepting gap in the 
circle traffic stream after checking the occupancy of some 
predetermined vehicle spaces around the circle (see Appendix 
FE). This crude approach of uniform gap acceptance procedure 
was chosen to avoid the prohibitive cost of gap and headway 
measurements in the field through the use of some more 
exacting procedures such as time-lapse photography. 

When the probability cf gap acceptance table has been 
constructed after analyzing field data from time-lapse 
photography, for instance, the table is employed in a gap 
analysis process by defining certain sections of the roadway 


which should be free of vehicles before a vehicle (which 
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happens to be evaluating gaps) would enter the circle 
roadway or the major street. This implies that the use of 
such procedures estimate the time gaps which are in turn 
translated into space gaps by drivers in the approach 
streams. Thus, by changing the gap acceptance logic in the 
model so that different drivers would define different 
safety zones, and hence, different numbers of vehicle spaces 
that must be free of vehicles before they would attempt 
circle entry, the use of the space gap acceptance procedure 


in the model would be more economical and easier to apply. 


5.3 Recommendations for Further Study 

The field validation study employed in the model is 
not that sophisticated, since most of the parameters were 
based on the gocd judgements cf drivers, thus eliminating as 
many traffic hazards as possible. More field validation 
study should be pursued to ccnfirm or more clearly establish 
the parameters that have been estimated in the model as a 
result of the crude field studies. 

Further study should be undertaken whereby the rest of 
the statistics generated by the model, such as the waiting 
times in approach queues, the mean system delay, and so on, 
would be validated by means of rigorous field studies. When 
most of the statistics as generated by this model, have been 
validated, it would then be possible to faithfully evaluate 


the suitability of the choice of traffic signal or traffic 
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circle as a method of intersection control for some specific 
intersection configurations. 

The model is built to accomodate a maximum cf 6 arms 
(or arteries) on a traffic circle, but the validation study 
was confined to traffic circles with 4 arms (the commonest 
form of a circle in the City of Edmonton). It will be of 
further benefit to the traffic engineer to have the model 
validated for 5- or 6-arm traffic circles so as to be able 
to evaluate and compare the performance of a S-way or 6-way 
intersection under a variety of intersection ccntrol 
methods. 

The number of lanes is set at 2 for all sections of 
the system. It will be worthwhile to alter the program 
Slightly so as tec be able to acommodate muiltilane 
specifications on different parts of the system 
configuration. 

Very often, the develcpers of individual intersection 
Simulaticn models mention extensions of their models to 
network simulation on intersection to intersection basis. 
One wonders how the traffic circles (which may happen to 
form part of the network of intersection) are treated in 
such network simulations. Now, it is possible to combine a 
refined form of the traffic circle simulation model and _ the 
existing individual traffic signal simulation models to 
realistically simulate a network of intersections on 


intersection to intersection basis. 
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APPENDIX A 


PROBAPILISTIC FUNCTIONS IN THE MODEL 


A.1 Unifcrm Distribution 

If integers are picked at random from the range A to B 
inclusive, the frokability of any particular one being 
picked is given by the density function 

f(y) = 1/B-A where A < y < B (1) 
The cumulative probability F is obtained by integrating the 
density function over the range of y 

y 
F(y) = Se cnvay 
a fa7B-a = y-A/B-A (2) 

and solving for y to obtain the relationship 

y = A+(E-A)F (3) 
From (2) we see that F varies between 0 and 1 Since y varies 
from A to B; so that for y in expression (3), we can 
substitute RN (random number from the range 0 to 1) and thus 
make y a random variable which is a function of RN: 

y = A+(B-A)RN (4) 
Denoting y by FN, we get the equation 

FN = A+ (B-A) RN (5) 
Equation (5) ccnverts values of RN to uniformly distributed 
values which lie in the range A to B as shown in Figure A.1. 
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Figure® Axl! 
Uniform probability distribution. 
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A.2 Exponential Distribution 


The exponential probatility density function is given 
by 

f(y) = exe (-dy)- (1) 
The cumulative probability is obtained by integrating (1) to 


give 


" 


y 
F(y) fexp (-hydy) 


1 - exp(-\ydy). (2) 


since y ranges from 0 toco, F varies from_0 to 1, so that if 
we denote a randcm number obtained from a U(0,1) as RN, then 
we can substitute RN for F and solve equation (2) to oktain 

y = 41n(1-RN) (3) 
The exponential distribution whose probability density 
function is that given by equation (1) has a mean of 3. If 
we denote this mean by m, and if we rename y as FN, equation 
(3) becomes 

FN = -m x 1n(1-RN) . (4) 
This is the cumulative exponential distribution function. 
Both the original density function and the cumulative 
distribution function are plotted in Figures A.2 and A.3. 
From equation (4), It 1s cleal that, the valucpor FN “can be 
obtained by finding the valve of 1n(1-RN) and multiplying by 


the mean m. 


o/- 
wae 


a 


tov te ai nodosns? yotedeb ysiiiiadang £61207 0qx0 


ty  ehh-panek = (yp 
a9 (t) peideapsrod gd bootetdo af qititdsdor; svitefosse 


a anni ins® > aes . war | 
: ™ res 
! 


odt 


eae = wt 

(S) -ieilaee “T= 
+s dedd oe <} oF 0 aoat coiTse a yeeod O eo7d eopnes ¥ enake 
eods eR ae {f.0)U & aoc? Bowistdo wedene aobast 6 etonsb 
aletis ot (S) acktsepe evioe Lan 4 tol @€2 ovusltedua asd 
i) (ta-tyale = ¥ 
ytiensb y2itedadcer4 ezodw netredéazely> leiseenogxe 
41 -+20 aaee 5s asd (1) ovtteueo yd aevite taedt aft nots 
fotteips \Wi.es y smeder Sv 22. bas 4s yi neon ald? esoneb 
aseooed 
(Pp) (Hi-fyal = o- = OS 
«tieiso au? wortudistete feitusesgqes oveatstiaeus off ai @ 
svitelowes odf foe achzoce? yrieaeh “Seateise ed? 
.t.4 fee S.4 eoawele af begtolq sas. avitooo. nottodt 
sd wed Wa k6 evlav gfe ted? festo ef +2 oth) noltsu pe 


yd potylqesion bas (¥8-F) al 50 eviev ott puibais yd bente 
. »@ nsoe 


a 


116 


Vispek Ve pice 
Figure A2 Exponential density function. 
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To use the exponential distribution to generate vehicle 
arrivals in a traffic model system, the INTERARRIVAL time of 
vehicles must be known. The interarrival time is the 
‘Inverse! of the frequency cf vehicle arrivals. 

For example, if 20 vehicles enter the system via a 
particular approach in a minute, then the frequency of 
vehicle arrival is 20 vehicles per minute and the 
interarrival time is 1/20 minute per vehicle or 3 seconds 
per vehicle. In cther words, one vehicle will arrive at 
this particular approach every 3 seconds. 

Without the specification of the exponential 

distribution, vehicles will enter the system at regular 
intervals of 3 seconds. 
The vehicles are made to enter the system according to a 
Poisson distribution by generating a number from the 
exponential distribution, EXRN , and multiplying by 3; in 
this way vehicle arrivals become random, and exponentially 
distributed. 

While describing the relationship between interarrival 
time and frequency, it should be pointed out that if a 
collection of interarrival times are exponentially 
distributed, the corresponding frequencies belong to a 
Poisson distribution. Hence the words "Exponential! and 
*Poisson' are often (but not always correctly) used 


interchangeably. 
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A.3 Normal Distribution 


A Normal distributicn is completely specified if its 
mean and standard deviation are given. 
The density function of a Normal distribution is denoted as: 

p(V) = 1/aro exp-[ (V -y) fo 2/2 (1) 
where 

V is a random variable 

Oo is the standard deviation and 

u is the mean. 

A standard normal distribution is defined as one with a mean 

of zero and standard deviation of 1. Its density function 

is obtained by substituting the values of 0 and 1 for and 
respectively, in equation (1) to obtain: 

p(F) = 1/427 exp(-F2/2). (2) 
Both the Normal probability curve and the standard form are 
Shown in Figures A.4 and A.5. 

A normally distributed variable V can be translated into a 
value F by using the definiticn of the standard normal 
deviate: 

ee uO. (3) 

===> V= OF tu. (4) 

If we can somehow obtain a value of F (belonging to the 
standard normal distribution) we can then use equation (4) 
to convert it to a value of V (belonging to the normal 


distribution defined by uw ando). 


ial Sh sn an ta ie 
van aodinived sasbuate Ban Anes 


ras botoneb al aokiudi 13028 en ae golsoau2 
("7 eNSE ae Wrase Om = oe 


oideiae? sobasx ‘s ok ¥ . 
Aes Be cost 
Bie dotestveh buahasse odd et dD 


¢ 

-asoa ods at u 

,ay Sbes 

econ 2 d¢iw eno at Bonkiab et aobsvdiareth {[sa100 = A 


avigosow? yibeae® 271 .T %o noiseived Gisbaste bae are - 
3 hee 
bes rol t bas 0 Yo esalae odd palsusizadus qd bentssdo a 


- 


)o as 

sntevdo ot {f) soiteepe ot thet oe ee 

(2) (8\58=)gxe FbN = (04 7 
» hoa Pet 


ote 702 brebsete of3 Doe SYIUD yrilidedory 1e2108 od? dtoe 
4 =i t@ 


-2.A Boo 8.4 nesopld at ors 

« offt bogsLens2? od ano ¥ eldetusy bopodiayals titeazoa - 

SES 

fewzon b$cebseta of? Yo apitiatieh si? ates or q oolsy 
; ‘Fy 


(€) oNurt y q 


(©) v4 > > = 0 Come 
' oe 

ade of pndonoted) ¥ 30 outer #, néssdo vodenoa B82 eM 

i ' a 


(4) netzoupe eeu nedt ma 9M iahevessteud Sanson baum 
Sauson sda ot pakpubledy ¥ 0 sofee 6 of #2 setter Od 
. (> be sxd Seatteb eoddodiasels 


Figure A 4 
Normal probability curve. 
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The cumulative standard normal curve is shown in two 
different forms in figures A.6 and A.7. The cumulative 
curve as shown in figure A.6 is a conversion from the 
Standard normal curve which yielded a monotone increasing 
function. From the curve it is seen that the probabilities 
associated with areas under the curve are cumulative so that 
there is a cumulative probability associated with every 
value of F, and that probability is equal to the area tc the 
left of the ordinate ccrresfending to F. 

By reversing the dependency of the variables and using RN 
instead of P, we get the curve as shown in figure A.7 The F 
will assume values ranging from approximately -4 to +4 as RN 
varies from 0 to 1, 

The use of the Normal distribution in discrete 
Simulaticn systems can be described as follows: the user 
wants to obtain random values V from a distribution with 
mean p and standard deviation o. To do this, he defines a 
cumulative standard normal distribution from which randon 
values F are obtained. These values of F are translated 
into values of V by means of a Simple equation as shown in 
(4). 

In traffic flow simulation studies, vehicle desired or 
target free flow velocity is usually determined by means of 


anormal distribution. 
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The mean for the velocities in a particular approach 
direction of flow is specified, and the individual desired 
velocities are generated normally about this mean. Usually 
the speed about this specified mean is limited within a 
specific range, so that some form of transformation has to 
be carried out if the generated random deviate falls outside 
the specified range. 

For example, if an arterial mean velocity is denoted 
by M, and target free flow velocities are generated normally 
about this mean velocity in the range of 0.75 to 1.25, then 
we have the inequality: 

OC ONeSmo VEL eS el 2 OM. 
In the inequality expression above TVEL is the computed free 
flow or target velocity for a generated or newly arrived 
vehicle at the specified generating point; so that for a 
mean arterial velccity of 24 miles per hour, the individual 
target free flow velocities will be selected from the range 
of 18 to 30 miles per hour. 

To obtain normally distributed velocities in the above 
range we require random valves V from a distribution of mean 
1,0emand “standard deviaticn “1/20. To do this, a standard 
normal distribution is defined with mean O and standard 
deviation 1 in the range -5 to +5. Then by the use of 
U(0,1) random numbers values of F (which are normally 
distributed with uU=0O and O =1) are obtained and translate 


them intc values of V belonging to the ‘nonstandard! normal 
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(5) 
procedure values of F are 
and use eguation (5) to 


Opvalmiivalnuessotevern the range 0.7 55to 1.25) and multiply eby 


M (the arterial mean 


arrival 


velocity) to obtain the 


individual target free flow velocities. 
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APPENDIX B 


COMPONENTS OF THE MODEL 


B.1 Vehicle-Index Arrays 

The physical system of the model is made up of three 
distinct sections namely: the approach, the circle and the 
exit sections. These sections are represented in the 
computer as follows. 

The apprceach and the circle sections are each made up 
of a two-dimensional array, the first of which is an index 
for the identification of a particular vehicle and the 
second dimension is an identification of the particular lane 
being occupied on the roadway. The exit section is made up 
of a one-dimensional array which stores the vehicle indices 
before the vehicles exit from the _ systen. The arrays 
themselves contain the storage addresses of the 
characteristics describing the vehicles occupying the two- 
dimensional roadway. 

As vehicles are intrcduced into the system, they are 
numbered sequentially and stored in the system in the 
approach array. During the simulation process, the movement 
of any particular vehicle is represented by updating the 
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position and velocity of the vehicle within the same column 
of. the matrix, or by shifting its index within the same 
matrix (when lane changing takes place within the same 
section of the system) cr by shifting its index from one 
matrix to another (as the vehicle leaves the approach and 
enters the circle or from the circle to the exit lanes) and 


the general register is updated accordingly. 


B.1.1 Approach Array 

The approach array * APPRCH({I,J) ' is dimensioned to 
accomodate a maximum of 50 queued vehicles, so that an 
intolerable situation will cccur in the system only when 
there are more than 50 queved vehicles in an approach lane 
at one instant during the simulation run. 
Vehicle processing in the approach lanes is done in a 
circular manner, so that the arrays are capable of accepting 
infinite number of vehicles so long as there are nc more 
than 50 vehicles in a lane at a particular scan interval. 
This circular array processing concept is described in 
detail in section B.2. 

At any time, an element of the approach array either 
contains a zero or a vehicle index aj 

APPRCH (I,J) = 0 if no vehicle occupies that 
section of the approach lane 
or = aj if a vehicle occupies that 


section of the approach lane. 
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The vehicle whose index is a; occupies the i-th position of 
the approach. lane j, but it is not necessarily the i-th 
vehicle in apprceach lane j. The i-th vehicle in approach 
lane j is the vehicle which occupies the (NFSTA(J) +i) -th 


position of approach lane j as explained in section B.2. 


Boi.2 Circle Array 

Thesmcitrcles, array §CIRCLEW 1 jd tsS dimensioned to 
accommodate a maximum of 10 vehicles per lane of a circle 
lane section. Unlike the approach lanes, where vehicles are 
always free to enter the lanes, vehicles from the approach 
lanes can enter the circle cnly when there is enough space 
in the appropriate circle lane (determined from some simple 
deterministic rules) so that lccking up of the circle never 
occurs. 

Like the approach array, processing of vehicles in the 
array is done by reference tc pcinters which point tc the 
first and last vehicles in the particular lane of the circle 
arLay. 

The odd-numbered lanes form a continuous inner circle 
lane and the even-numbered lanes form a continuous outer 
circle lane. 

Thus, a vehicle in the inner circle lane travels fron 
one inner circle lane secticn to the other until it comes to 


the point of, exiting from. thescircle. 
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Boles UX ite Array 

The exit lanes are represented by a one-dimensional 
PEAY sok os eee This is because exit flows in the system 
are recognized only as one vehicle space, and that a vehicle 
is removed from the system after travelling one vehicle 
space along the exit lane. 

Vehicles from inner circle lanes exit from the system 
via the inner exit lanes, and vehicles from the outer circle 
lanes exit from the system via the outer exit lanes. 

The system is organized in such a way that vehicles in 
the circle lane J always exit from the system via exit lane 
J; so that at any moment in the simulation, the element e, 
of the exit vector contains either a zero or the index of 


the vehicle occurying the exit lane J. 


B.2 Vehicle Lists 

In accordance with the circular list processing 
concept in the model of Gerlcugh and Wagner [19], vehicle 
lists are maintained in this model in the approach and 
Gircile: .SecuLons sot the sisysten. Unlike the orthogonal 
intersection of Gerlough's model where an approach lane 
Goulds forms ae continucucmlis twi the anbwiexat®, Gane je in erthis 
model, approach lists are taintained separately from circle 
Lists. 

Lists are not maintained on the exit lanes, since 


vehicles are deemed to have left the system aS soon as they 
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cross one vehicle length (vehicle space) along the exit 
lane. Approach lists start from the system entry boundary 
and end at the circle entry line. Each approach lane in the 
system forms a separate list, so there are as many approach 
lists as there are approach lanes. Similarly, circle lists 
are maintained in each lane section of the circle depending 
on the number of lanes arcund the circle; so there are as 


many circle lists as there are circle lane sections. 


B.2.1 Approach Lists 

Vehicle lists on the approach are ordered so that 
vehicle sequence on a list is identical to vehicle sequence 
on the roadway. The first vehicle on a list has the highest 
position value (that is, farthest from the zero or reference 
point), and the last on the list has the lowest position 
value (closest to the zero cr reference point). 

When a vehicie is generated at an approach boundary 
(hereinafter called ‘System Boundary"), some vehicle 
parameters are attached to it, among these parameters is the 
position parameter. When the vehicle is generated it is 
placed at the bottom of the list by giving the value of zero 
to the pesition parameter; and the movement of the vehicle 
on the approach is done by incrementing the value of the 
position parameter until it has travelled the full approach 
distance. 


Thus, when a vehicle is generated and placed cn the 
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approach array, itsS position is not physically shifted 
around until it is ready to merge with the circle traffic 
stream. AS a vehicle merges with the circle traffic stream 
it is struck off the top of the approach list, and placed at 
the bottom of the anpropriate circle list depending upen the 
Originating approach lane which was occupied by the merging 
vehicle. 

Control over vehicle processing on a given approach 


lane is maintained by reference to the following index 


arrays: 
NFSTA(I) = Index of the first vehicle 
On approach list I; 
NLSTA(I) = Index of last vehicle 
on approach list I; 
NA (I) = Number of vehicles 
on approach lane I. 
Thepindices of theyfirst and Jast vyohiches on) ethe!l,.arproach 


lists are kept in separate arrays as indicated above, so 
that during updating, scanning of the list could be limited 
only to indices within the range of NFSTA(I) to NLSTA(I) . 

Initially, before any vehicle has been generated into the 


system the index arrays are initialized as follows: 


NFSTA(I) <----- ts 
NLSTA(I) <----- 0; 
NA(I) <--------- 0. 


As the first vehicle enters approach lane I, it is 
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identified as the first on the approach by giving it a 
number equal te 1 on the list; and by storing the index of 
its position on the approach list in the vector NFSTA(I) to 


give the following status of the arrays: 


NA(I) <-------- NA(I) + 1 = 12 
After the tenth vehicie has entered the approach lane I, the 


indices cf the approach list attain the status: 


NFSTA(I) = 1; 


NLSTA(I} = 10; 


NA (I) 10. 
We see that until the first vehicle on the approach lane [I 
is able to merge with the circle traffic stream, the value 
of the index NFSTA(I) stays at 1, irrespective of how many 
vehicles have entered the system via that approach lane. 
However, that of NLSTA(I) keeps increasing as more and 
more vehicles enter the system via approach lane I. As the 
Pirstevent€iSlassakle> CSPPaangeW**eiekt pts) circle *raffic 
stream we get the following Situation of the vehicle index 


array for approach lane I: 


NFSTA(I) <----- NFSTA(I) + 1 = 2; 
NLSTA(I) <--------------- 10 8= 410; 
NA(I) <-------- NA(I) - 1 = 9. 


The second vehicle cn approach lane I becomes the first on 
the approach as. the first vehicle on the list leaves; and 


the number of vehicles on the list is decreased from 10 to 
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9; however the tenth vehicle on the list still remains as 
the last on the list until another vehicle arrives on that 
approach lane. 

In summary, as more vehicles enter approach lane I, 
the value of the index NLSTA(I) is always incremented in 
order to point to the most recent vehicle on the appfroac? 
and the number of vehicles on the approach list is 
incremented accordingly: 

NLSTA(I) <----- NLSTA(I) + 1; 

NA(I) <--------- NA(I) + 1. 

Overflow of vehicles occurs on approach list I when either 
of the fcllowing conditions cccurs: 
NLSTA(I) > JAPP and NFSTA(I) = 1 
or NLSTA(I) = NFSTA(I) and NA(I) = JAPP 
where 
JAPP is the maximum number of vehicles that can be 
accomodated on an approach lane conveniently. 
However, if NLSTA (I) > JAPP and NFSTA(I) # 1 
then NLSTA(I) <----- 1. 
Similarly, as more vehicles leave the approach lane I, the 
index vectors are updated as follows: 

RESTA (LW <sasts | NESTAr(L OA tiel 5 

NA(I) <-------- NA(I) - 1. 

When NFSTA(I) > JAPP, NFSTA(I) is automatically reset to 1 
whether or not there are any more vehicles on the apfroach 


lane. The number of vehicles on the approach list I, at any 
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instant, is given by the value of NA({I) which is ufdated 
continuously during the simulation. However the value of 
NA(I) could be computed at any time from the formulas: 

NA (I) = NLSTA(I) - NFSTA(I) + 1 if NLSTA(I) > NFSTA(I); 

or = ( JAPP + NLSTA(I) ) - NFSTA(I)+ 1 
if NLSTA(I) < NFSTA(I) . 
The pointers NLSTA and NFSTA are reset to their initial 
values whenever there is no vehicle on the appropriate 
approach list. | 
That is, if NA(I) < 0 then NFSTA(I) <----- 1 
and NLSTA(I) <----- 0. 

Thus, the pointers NFSTA(I) are moved(incremented or reset 
to 1) when a vehicle leaves an approach lane and merges with 
the circle traffic stream; and the pointers NLSTA(I) are 
moved (incremented or reset to 0 or 1) as more vehicles join 


the appreach lane I. 


Be2.2eCircle Lists 

Vehicle lists on the circle lanes are treated much the 
same way as vehicle lists cn the approach lanes. Index 
arrays NLSTC(I), NFSTC(I) and NC(I) are maintained to 
monitor the movements of vehicles around the circle in the 
same way as vehicle movements on the approach lanes are 
monitored by references to the index arrays 
NLSTA(I) ,NFSTA(I) and NA(I) as described in section B.2.1. 


But there is more shifting of vehicle indices around 


‘ (ayacene > (ayarege te oe: owes 
{siting stiedt of IeRB2  9T6 Yat bas a€aa-exesmtog wd? . 
etwizgeagqs of% no slokdee of ef erect 1oveasdy espisv 
: * | «tet dosozggs 

po mas-- > (Spatew aede @ > (TpAW Bt yal sad? 

0 -----> (ipo Bas pwd tewed 

tetsa to haeimenetont) bevo8 ems (tY ATEN ex9satog eds Rude 
d¢tw ssostee bos onset fosdagge ne sovaal ofotiev 4 aedv (f oF 
ots (L)ATSda azerhiog of bnew iwtedte sbitext oforio ‘edt 
ntot seloidew s20n as (f 20 0 oF +ee67 29 bednemezcat) Sevos 
«T sask doso3qqs’ ode 

see 

ase21 olpal> SiS 

add dowe bosaunt ers asian elorto edt do-eteni elobiev ~ 
sebat = 6-@ensi dpsongge sit a> edebl Siotiaw es yaw ease 
o¢ baaisinise e175 (T)98 Sas (EpOTETM. 4 (1)pOTeMH) (eyeure 
ait gt ofosko os, bavoas) cefobisl 46 e¢nembyon oft) to¢zaom 
ean axanl dosoaqys. sis no 6 etnem@syot siotiey as ys Sasa 
eyes «awh © oto DUSSTH2 «gk sasOsAMON 
sFis@ noltos® ob Gaditzsued ab (BV AM bus @yavenn, (xparege 
bauors apotbat Lottew Io gabdhkde stom ek ssede 908 | 4)" ! 


—— 


138 


the circle than at the approach lanes. Unlike the approach 
flows where vehicle physical position on the approach array 
stays the same until it merges with the circle; vehicle 
index on the circle is shifted from one vehicle list to 
another as the vehicle completes its journey in one lane 
section and travels to the next lane section, until it 


finally joins the exit lanes. 


B.3 Vehicle-Characteristic Array 

The characteristics of all vehicles in the system are 
stored in each row of a vehicle characteristic array, 
AUTO (I,J), so that pertinent information about a vehicle is 
obtained by reference to this array. When a new vehicle is 
generated at an approach boundary (also called the Systen 
Boundary), its characteristics are stored in a row cf the 
characteristic array and the row number of this array is 
stored in the appropriate approach lane. During simulation, 
it is this rew number (rather than the vehicle sequence 
number) which is shifted arcund in the system to represent 
movement of the vehicle. 

Whenever an updating routine is to be performed ona 
particular vehicle, its sequence number on the vehicle 
characteristic array is extracted from one of the vehicle 
index arrays; and the updating process is carried out by 
modifying the appropriate elements of the row cf the Auto 


array corresponding to that number. In the beginning cf the 


as co i. i 
ga 
eividey yelotto ed? dthy espree 9f Lisun ones) of: ay . 
se kta in nie alae ale ie ela a 
anel ono of yoatwot ast eorotgats. ofottey au? 2B zedso05 . 
+4 Lisen wottose aust faen e4% 09) ahever? bas” sosisen < 
.20niL Se 063 cote: Yisats 
met vesngey, . 


- 1 © 
(ers SPgetretosrsds-pfotdev Ess 


ers eeteye ott at eslotdew Dis to sobteitesosisdo odt?- asig ; 
,{2its aiveliezoeqsd> elotdew 6 I) vos dose at bezota 7 
si ef{vider & teods aotrnezc30 + neakdreq dads OF v (ST) O8SE | 7 
ai sioiftev eon s aed .warto Sida 6% 90nsTa2e2 yd beatetdo — 
sevay2 soo Seliss onks) ¢zehbavod tosorggs a5 26 bosszeneR 7 
ait 29 403 © a¢ heicta eta sobvelretostsda ett rer i 


el yarxs atdt Yo wedudia wor odd bos ye7rs otvetretoersis 


’ 
,solssinaie paized .enal dasoig7s etekagoaqys eds at beset 7 - 


scaswpee einidev of? asd? Jedtet) tedeon #7 sidt at. yt 
a of eovoys wit at Banos eaRbte aL dtde tocune 
Levniesprimatiaat > 
» #0 bosrcdaeq st of et snizvoy pehsenqy as mennaneth yAROy 
cividey sft 40 1adNon sbasupee ert yess toot s3q 
einidey ad? 20 .9H0 sort bemmeaiem pe yeTse. 


_ 


- ee| 


134 


Simulaticn, all rows of the vehicle characteristic array are 
linked together, and an index of the first row is stored in 
the link pointer, LINK. When the first vehicle arrives in 
the system its characteristics are stored in the first row 
of the vehicle characteristic array. Thereafter, 
characteristics of subsequent vehicles are stored 
sequentially in the rows cf the array until some cf the 
previously occupied rows are vacated by exiting vehicles. 

When a new vehicle arrives in the system and a 
previously occupied row of the vehicle characteristic array 
is vacated, then the link pcinter would be pointing tc the 
most recently vacated row, which will then be assigned to 
the newly arrived vehicle. Thus, arriving vehicles are 
always assigned to a fresh Space in the vehicle 
characteristic array only when there are no vacated spaces. 
This form of ‘Garbage Collection procedure in using the 
vehicle characteristic array insures an optimum use of the 
array; and in this way , saturation of the system (when 
there are at least 1000 vehicles in the system) will occur 
only when the link pointer points to the last row cf the 
vehicle characteristic array. 

Typical parameters stored in the vehicle 
characteristic array for each vehicle which enters the 


system are shown in figure B.1. 
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The numbers in Figure B.1 indicate the columns of a typical 


re ee ee | eee 


Figure B.1 Vehicle Parameter List 


row of the AUTO array. 


4 


10 


11 


vehicle sequence number, referred to 

as vehicle identification number;. 

entry lane of the vehicle, which of necessity 
must refer tc an arproach lane; 

exit artery of the vehicle; 

time of entry into the system; 

the total distance to be travelled 
around the circle; 

the speed of the vehicle 

at any time during the simulation; 
position of the vehicle in the system 

at any time during the simulation; 
vehicle desired velccity; 

time of joining a gueue at the approach; 
time of entry into the circle; 


time of exit from the system. 
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APPENDIX C 


PROGRAM LOGIC AND INPUT ORGANIZATION 


C.1 PROGRAM LOGIC 

The logic of the fprogram is explained in detail in 
this section with accompanying flow charts. 
The order of the logic segments in this section does not 
necessarily duplicate the order of processing in the model 
as explained in Chapter III. The following abbreviations 
have been made use of in the accompanying flow charts. 

Q-IN = Inner approach lane queue; 


Q-OUT 


Outer approach lane queue; 


C-IN-R = Inner circle lane section to the 


right of a potential merger; 
C-OUT-R = Outer circle lane section to the 

right of a potential merger; 
C-IN-L = Inner circle lane section to the 

left of a potential merger; 
C-OUT-L = Outer circle lane section to the 


left of a potential merger; 
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IN-ZONE = Conflict zone of a potential merger ready 
to merge from Q-IN; 
OUT-ZONE = Conflict zone of a potential merger ready 
to merge frog Q-OUT; 


RF 


Random fraction generated from U(0,1); 


AFRQ Arrival frequency at an approach direction; 
Potential merger = Leader of approach gqueue(or lane) ready 
to accept lag/gap to enter circle. 


Some of the abbreviations above are illustrated in Figure 


C.1.1 Vehicle Generation 
The assignment of origin-destination is controlled by 

the user, who sets the allowable origin-destination 
distribution at an input boundary. The flow rate, indexed 
by boundary is also provided by the user. 
If the volume of arriving vehicles is denoted by VOL 
vehicles per hour, then in an interval of time equal to the 
scan interval, the average frequency of arrival, denoted by 
AFRQ is given by: 

AFRQ = VOL x SCAN/3600 
where 

SCAN stands for the unit of scan interval; so that 
for a scan interval of 1 second, the average frequency of 
arrival is given by: 


AFRQ = VOL x 1/3600. 
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A random fracticn with rectangular distribution about a mean 
of 0.5 (that is, U(0,1) random number) is generated for each 
scan interval and if the random fraction RF is less than or 
equal to the average arrival frequency AFRQ , then a vehicle 
is assumed to have arrived at the particular vehicle 
generating point during that scan interval. 
However, if RF > AFRQ then no vehicle is assumed to have 
arrived. 

Vehicle arrivals then follow a Binomial distribution 
with a constant probability which is an approximation to a 


Poisson distribution. 


C.1.2 Assignment of Directicn 

If a vehicle arrives in the scan interval, then it is 
assigned a directicn by the Monte Carlo technigue. A random 
fraction, RF, is generated and compared with the estimated 
probability of a vehicle turning right, turning left or 
going straight through the intersection. 

The probability of a right turn (RT) is taken as being 
equal to the percentage of right-turning vehicles over a 
period of time, and if RF < RT the vehicle is assigned a 
right-turn destinaticn. 

Similarly, the probability of a left-turn (LT) is 
taken to be equal to the percentage of left-turning vehicles 
over a. peniod...of ,.time,,. and, if. RF .\s. LT, the, vehicle is 


assigned a left-turn destination. 
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Finally, the probability cf a straight-through (ST) is 
taken to be equal to the percentage of straight through 
vehicles over a period of time, and if RF < ST the vehicle 
is assigned a straight through destination. 


The probakilities RI, LI and ST are cumulative in the 
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C.1.3 Assignment of Lane 
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the inner approach lane as is denoted as ‘lane 1 and the 
cuter approach lane as ‘lane 2', then all right-turning 
vehicles are assigned to the right-turn queue (lane 2) and 
all left-turning vehicles are assigned to the left-turn 
queue (lane 1). 

Of the straight through vehicles, assignment to lane 1 
or lane 2 is accomplished in cne of two methods as explained 
below: 

1. The lengths of the two queues of the approach lanes 


are compared and the straight through vehicle is 
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assigned to the shorter queue. 

2. If neither of the queues appears shorter than the 
other, then the assignment of the vehicle to lane 1 or 
lane 2 is predicted with the Monte Carlo technique as 
follows: 

A random fraction, RF, is generated and compared with 
the probability based on observed lane usage of 
straight through vehicles from the approach direction, 
and the vehicle is assigned to the correct lane in 
that random fashion. 
Presently, it is assumed that vehicle usage of inner 
approach lanes by straight through vehicles from all 
approach directicns cf the traffic circle is 0.70 (p) and of 
the outer approach lane is 0.30 (1-p). Of course, since p 
a is read into the program aS an input variable, the 
probability of lane usage by straight through vehicles from 
the various afrfroach directions can be set to any desired 
value based on the data at hand. See Figure C.2 for the 


flow logics of sections C.1.1, C.1.2, and C.1.3. 


cC.1.4 Circle Aprroach Flow 

The logic for processing vehicles which are travelling 
from the system boundary to the circle is the same for both 
inner and outer approach lanes. The difference comes only 


in the logic for gap acceptance at the circle. 
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Essentially the logic fcr processing vehicles on the 
approach lanes (illustrated in Figure C.3) as follows: 

The leader of the approach lane is first considered. If the 
leader is within some specified distance (presently taken to 
be 15 feet) from the circle, then the gap acceptance logic 
is called to determine the probability of the vehicle 
accepting the available gap in the circle traffic stream. 
If it accepts the gap to enter the circle, it is processed 
out of the apprcach lane into the appropriate circle lane, 
and the necessary statistics are gathered to update the 
status of the originating approach lane and destination 
circle lane. 

However, if the leader rejects the available gap or if 
there is no gap in circle stream, then the vehicle is 
decelerated to a stop at the circle entry line within the 
next scan interval. If the leader is further away than the 
specified distance (here, 15 feet) from the circle, then it 
is accelerated or decelerated depending on how far it is 
from the circle and its current speed. But the leader is 
not allowed into the circle during the current scan interval 
if it is further away from the circle than the distance as 


specified above. 
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After the leader has been updated, the rest of the vehicles 
in the approach lane are updated by merely computing the 
distance headways between successive vehicles and assigning 
new velocities £oxr the next scan interval. After 
reassigning the velocities, the position parameters of the 
vehicles are updated to show the respective positions of the 
vehicles during the coming scan interval. When all vehicles 
of a lane have been updated as described, the vehicles in 
the next lane of the current set of lanes are considered in 


the order as indicated in Chapter III. 


C.1.5 Gap Acceptance from Inner Approach Lane 

For a vehicle wishing to merge with the circle traffic 

stream from an inner approach lane, it is faced with the 
evaluation of lags and/or gaps from the two circle lanes 
(circle eae ee Since its path to its destination in the 
circle crosses both circle streams of traffic. 
Such a vehicle is concerned with that half of the traffic 
circle which immediately faces it, and as far as its 
evaluations are concerned, the cther half of the circle does 
not exist. In its immediate half of the circle there are 
two continuous circle lanes (inner and outer). 

To distinguish between destination lanes and gap 
evaluation lanes in the circle, the two continuous streams 
facing the potential merger are broken up into two halves 
each. (see Figure C.1). One set of lanes is from the 


entrance of the approach artery on his immediate left to his 
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own exit artery (C-IN-L and C-OUT-L), and the other halves 
are from his artery to the exit leg of the artery on his 
immediate right (C-IN-R and cC-OUT-R). To the potential 
merger, we can therefore distinguish between right and left 
circle lane secticns. His conflict zone happens to be the 
first vehicle pcsition of the outer and inner circle lanes 
to his right. 

Before such a vehicle attempts an entry into the 
circle, the fclleowing conditions must be tested and 
satisfied simultaneously: 

(1) that the last vehicle in C-IN-R is beyond the conflict 
zone; 

(ii) that the last vehicle in C-OUT-R is beyond the 
conflict zone; 

(iii) that if the first vehicle in C-OUT-L is not an 
immediate turner (that is, not turning at the exit of 
the potential merger's artery), then such a C-OUT-L 
vehicle should not reach the conflict zone in the 
outer circle during the next scan interval; 

(AV\nthat Vit the first |) vehicie in’ C-1IN-L) is» not) an 
immediate turner (that is, not turning at the exit of 
the potential merger's artery), then such a C-IN-L 
vehicle should not reach the conflict zone during the 
next scan interval. 

When the four ccnditions abcve have been met, the potential 


merger is said to have accefted the gap at the circle. 
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C.1.6 Gap Acceptance from Outer Approach Lane 

The gap acceptance logic for vehicles wishing to enter 
the circle from the outer approach lanes (Q-OUT vehicles) 
differs from that of the potential mergers from the inner 
approach lanes (Q-IN vehicles) in the sense that, the former 
evaluate gaps in the outer circle lane sections only. 

A Q-OUT vehicle considers C-OUT-L and C-OUT-R lanes in 
much the same way as a Q-IN vehicle considers the same set 
of the lanes; and accepts lag and/or gaps to enter the 
circle when the necessary tests (incorporated in Section 
C.1.5) are satisfied. 

One important gap acceptance feature for Q-OUT 
vehicles is the fact that they can froceed simultaneously 
into the circle as soon aS a Q-IN vehicle accepts a gap and 
enters the circle. But the reverse is not true. This is 
because the gap acceptance logic for Q-IN vehicles 
incorporates the gap acceptance logic for Q-OUT vehicles but 


not vice versa. 


C.1.7 Inner Circle Flow 

Within the system, vehicles in the inner circle lanes 
have the highest pricrity whenever there is any conflict in 
movements. The general flow logic in both circle lanes is 
the same as that in the approach lanes; in the sense that 
during a scan interval one and only one vehicle (the leader) 


is processed cut of the circle lane list into the next 
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circle lane list or into an exit lane (see Figure C.4). 

The rest of the vehicles in the lane list are merely updated 
in the same manner as the rest of the vehicles in the 
approach lanes. 

The leader of an inner circle lane section either 
travels to the next inner circle iane section or exits from 
the system when it is within some specified distance (15 
feet) from the end of the current lane section. During 
exiting, an inner circle vehicle crosses the path of an 
outer circle vehicle with absolute priority. Since lane 
changing within the vicinity of the circle is not allowed in 
the model, it is during an entry into or exiting from the 
circle only, that an inner circle vehicle ever crosses the 


outer circle lane. 


Cs1.8 Outer Circle Flow 

Vehicle flow in the cuter circle lanes is the same as 
that for the inner circle lanes except that vehicles in the 
outer circle lanes have to yield right-of-way to vehicles in 
the inner circle lanes when conflict in movements occurs. 
After entering the circle from an outer approach lane, a 
vehicle in the outer circle lane travels from one outer 
circle lane section to the other until it reaches the 


exiting artery where it drives into the outer exit lane. 
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Figure C.4 
Inner circle flow. 
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During its passage from its current lane section to the next 
(anticlockwise), it checks fcr possible conflict with an 
adjacent inner circle lane vehicle which may be ready to 
exit. If such a conflict exists, then the outer circle 
vehicle yields right-of-way to the inner circle vehicle. 
When all vehicles of an outer lane section have been 
updated, vehicles in the next outer circle lanes are 
considered until all cuter circle lane sections have been 


updated (see Figure C.5). 


C.1.9 Exit Process 

There is no list processing along the exit lanes, 
because vehicles are allowed to enter the exit roads from 
inner and outer circle lanes one at a time, and moreover, 
vehicles are deemed to have left the system after traversing 
ane vehicle space on the exit road. After travelling one 
vehicle-space on the exit road, the following statistics are 
gathered about the outgoing vehicle (if the warm-up feriod 
for the run is cver): 

1. The time of leaving the system; 

2. Originating approach lane; 

3. Exiting artery; 

4, Desired travel time; 

5. Actual travel time; 

Ge pelay: 


7. Exiting speed. 
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When the above statistics have been collected, the vehicle 
is removed from the system Ly discarding its sequence number 


on the vehicle-characteristic array (see Figure C.6). 


C.1.10 Queue Measurements 
Queue lengths are updated for each approach queue 
during each scanning cycle. Initially, because there are no 
vehicles in the system, queve lengths are set to zere for 
all n approach lanes ( KUELEN(I) = 0, I=1, n), and the 
updating routine is accomplished as follows: 
(a) If KUELEN(I) = 0 and the first vehicle of this approach 
lane I (whose index is given by J = NFSTA(I) ) stops, then 
KUELEN (I) <----- A 
(b) If KUELEN(I) > 0 then 
(i) We decrease KUELEN(I) when 
the first vehicle on the approach 
which is not yet a member of the queue 
(whose index is given by J=NFSTA+KUELEN (1) ) 
Stops: KUELENAS) i<ce--- KUELEN (I) +1. 
(ii) We decrease KUELEN(I) when the 
first vehicle in the approach lane 
(whose index is given by NFSTA(I)) is 
able to merge with the circle traffic 


streams KUELEN(I) <----- KUELEN (I) -1. 


ey a.) 
wy, 

cet : aa 
; 10S SSA | 
ginidey off ,beddeLivo awed 


(3.9 sauptt soa) 


a ee 
ps ie 


“ohee _ 

=| 

 eonid 
edt A 


7 | | : ; Ae or.t 7 7 me : 7 
ee | ie me 
eimbp so80Rqgn dons 03 lilies 8 as adeno sae “4 
on e168 sus? pevsoed ,yliessiat ott sega e 0b 


102 oes oF J8a ate adtenel avong : eat 
we bas y(n ,f=t .0 = Gyuenson  ) wana soebagys 
renoifo aa bedeniqaggas of onktuor 
donviqqa eid? Yo sieldew a6ze9 ond pee O'S ywaaaon aI 
ned? yeqota ( (Tpaseaw = bb yw vex et sebad sod) I 
Hoseao> (0) WEIETA 
weve 0 < (TpMeIEUA BT 4 
‘adv (TYMRIRO® seeeI9e ON (2) 

ipserqga ois ao stotdew dats oat 

sueag edt To nedeem 5 soy oka 


(TyMaiawarATeTMEG Yd oew hp a rebut - 


fe (Lyuataue 


anal donorggs oad: me 
at (CR) AaeeH yd aBvER ak HOBML geome 
vantess vabonke tt am opzan “et ofan 


o xi a a iy ay i" i? 


Los 


Consider an | 
exit lane. 


Compute and gather the 
following statistics: 

!. time leaving system. 
2.0originating opproach 
.exiting ortery. 
4.desired travel time. 
.actual travel time. 

6. delay. 

7 total na of cars for this 
origin-destin. route. 


Remove car 
from exit lane 
& from the 


Figure C.6 
Exit processing 


154 


Using this technique of queue measurement (as illustrated in 
Figure C.7), any vehicle that joins a queue by stopping 
behind the last vehicle in the queue continues to be a 
member of that queue until it merges with the circle traffic 
stream. Thus, at any instant, all vehicles counted as 
members of a queue are not necessarily stopped; and they 
continue to be members of the queue from the instant they 
stopped behind the queued vehicles to the instant they are 


able to merge with the circle traffic stream. 


C.2 INPUT ORGANIZATION 

To be useful as a tcol for the traffic engineer in his 
evaluation study of a traffic circle performance, the 
acquisition of data to be used as input to the simulator 
should not constitute another problem in itself. 

Taking into account the types of data on the traffic 
circle that are generally available to the engineer, and the 
types that are readily collected if not available, there are 
3 different ways in which data would be accepted by the 
simulator depending on the type of manipulations that the 
engineer wishes to do with the available data. The 
following general form of data is usually available to the 
engineer; since data collected by field surveys on the 
traffic circle is readily and easily analyzed, using the 
traffic circle data analysis technique as outlined in 


Appendix D. 
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(a) Total volume (in pericds of 15 minutes, 30 minutes or 
1 hour) inflew of vehicles from each approach 
direction to the traffic circle. 
(b) The turning volumes from each approach directicn for 
the period of data collection. 
These two sets of data are always read into the program, so 
that they could be varied as often as needed. 
Now, the rest of the input phase depends on what the 
engineer would want to do with the above general form of 
teat tyes scireclerm@data® As already mentioned, the simulator 


accepts all three types of data format. 


C.2.1 Data Type 1 
In this type, estimates are made of the percentages of 
the total volume of vehicles that arrive at the approach 
directions at discrete time periods of t minutes 
(t=5,10,15,20,30 or even 60 minutes) assuming that the 
turning rates stay the same over the Simulation runtime. 
Therefore the time periods and/or the arrival 
percentages at the approach directions could be varied, and 
hence vary the arrival rates at the approaches during the 
specified simulation runtime. 
This form of data manipulation for different sets of time 
periods and different sets cf arrival rates can be rerfeated 
as often as needed. 
For example: 


Suppose that the arrival rate at the approach 
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directions of a traffic circle is varied every 15 
minutes fcr the first simulation period of 1 hour; and 
for the next hour of simulation the arrival rate is 
varied every 20 minutes. 
The number of times that the data is to be manipulated 
in this manner will have to be specified and read in 
as input to the simulator. 
Consider the data for the East approach as shown in Tables 
Getanancsaz. 
The hourly volume of arrival = 865 vehicles per hour 
therefore, on the average the frequency of arrival is given: 


AFRO = 86573600 


0.24 
Considering the percentage arrivals during the four 
15 minute periods of the simulation we get the following 


distribution of discrete arrival rates: 


Time Periods Arr. & Volume Avg. Arr. Freg 
430--445 25 216 216/900 = 0.25 
445--500 31 268 268/900 = 0.30 
500--515 24 207 207/900 = 0.24 


BOR a 20 173 1737900 = 0.21 
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By plotting the graph of arrival frequency against time 
periods, we find (in Figure C.8) that the periodic arrival 
distribution gives a better indication of the traffic 
fluctuations during peak hour flow than the constant average 


arrival frequency. 


Preparing input for Data Type 1 

The data format for Data Tyre 1 is shown in Tables C.1 and 
Ci2. In order to manipulate data with the simulator as 
outlined above the data should be arranged in the following 


sequences: 


(i) CARD 1 


should have the number '1' in column 4 


(ii) CARD 2 

should have the runtime (30 minutes or 1 hour) in 
seconds, punched in columns 1-4 right justified. A one hour 
simulation runtime will be punched in the first 4 columns as 


3600. 


(iii) CARD 3 

should have 1, 2, 3, ~«..-- Punched in columns 1-4 right 
justified, depending on the total number of sets of data to 
be manipulated in this fashion, or how many times the 
engineer wishes to manipulate the given data using this 


format. 
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Figure C. 8 
Arrival frequency distribution. 
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(iv) CARDS 4-6 


CARD 4 --- contains LFI volumes from all approaches; 

CARD 5 --- contains THR volumes from all approaches; 

CARD 6 --- contains RGT volumes from all approaches. 
(v) CARD 7 


contains the total arrival volume for the arrival 
arteries; and are punched in 6-column fields right- 


justified. 


(vi) CARD 8 

contains the time interval representing a period in the 
Simulation after which necessary generating characteristics 
are varied. This time must either be a multiple of 30 (5, 
10, 15) or a multiple of 60 (5, 10, 15, 20, 30, 60) or even 
a multiple of 120 depending on the specified simulation 


runtime on CARD 2. 


(vii) CARDS 9-? 
depending on the interval chosen in (vi) and the runtime 
chosen in (ii), there can be as many number of cards for 
this input phase as necessary. 
For example: 
(a) For a runtime of 1 hour = 3600 seconds anda 
15 minute period cf varying the generating 
characteristics we get the following 


INT1 = (3600/760)/15 = 4 


INT 2 INT1 + 1 = 5. 
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(b) For a runtime of 30 minutes = 1800 seconds anda 

15 ninute period cf varying the generating 

characteristics we get the following 

INT17 = (1800/760)/15 = 2 
INT2 = INT? + 1 = 3. 

In example (a) there will be 5 cards for inputting the 
percentages: CARDS 9-13. And for the example (b) there 
would be 3 cards for inputting the percentages: CARDS 9-11. 
In all cases, the first card, that is CARD 9, would contain 
the percentages of the total arrival volumes used for system 
fill-up purposes; and the rest of the cards should contain 
figures (integers) which add up to 100 columnwise. The 
cards are punched in 6-column fields right-justified. If 
there are more sets cof data cf this form in the sense of 
period to be used and the percentages to be used; then 
different sets of CARD 8 and CARD 9-? are repeated, but 
the number of such data should not exceed the figure as 


represented in CARD 3. 


C.2.2 Data Type 2 

In this type of data organization, estimates of the 
turning percentages from the approach directions are made at 
discrete time periods of t minutes (t=5, 10, 20, 30, 60) 
assuming that the average arrival frequency stays the same 
over the period of the simulation. 


The time periods and/or the turning percentages from 
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the approach directions are varied so as to have variable 
turning rates from the verious approach directions during 


the simulation period. 


The data format for Data Type 2 is shown in Tables C.3 and 
cC.4. The procedure for preparing input data to manirfulate 
the generating characteristics as outlined above is as 


follows: 


(i) CARD 1 


should have the number '2' in column 4 
(ii) to (iv) are the same as for Data type 1 


(vii) CARDS 9-? 

depending on the interval chosen in (vi) and the runtime 
chosen in (ii), there can be as many number of cards for 
this input phase as necessary. 
For example, if there are 6 periods for varying the turning 
percentages, then together with the warm-up period there 
would be 7 cards for each approach, the first of which is 
meant to be warm-up turn percentages for each approach 


direction. 
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i A A PES SUSE ee kk See aa 
TURN NOLUMES {TURN FREQUENCIES | 


{APPROACH | TOTAL - { 

{ { VOLUMES| LFT THRU RGT{ LFT THRU RGT | 

-————___+-__-_-___—__+-----_-_-_____—___-+1 me 

| 1 { 815 { 331 404 80{ 0.40 0.49 0.09 | 

{ { { | { 

{ 2 | 788 (9859 8657 46, 0.10 0.83 0.05 | 

{ | { { { 

| 3 | 988 { 265 680 G3/70526 0268 0208 { 

{ { ! { { 

| 4 { 865 ip Sree 70 801 90210°-0..89" 0.097 | 

(bpd re moe he J 

Table C.4 
Sample Data Format 2 
Periodic Turn Percentages 

Cele ti ee Ce a ae a Rear tre saat NETS AEE Te A Ea Te 
{15 MINUTE| SOUTH | NORTH | WEST { EAST { 
{| PERIOD { APPROACH { APEROACH {| APPROACH | APPROACH j 
{ {| (STURN) { (%TURN) { (ATURN) { (*TURN) | 
{ {LFT THR RGT|LFT THR RGT{[LFT THR RGT|LFT THR RGT{ 
asa. a laa eee See peste es eee ts 
( { | l { { 
{| 439-445 | 30 40 304 15 14 9tp149-67 9f'e10re305260] 
| { { | { { 
(R4*S—500 [e2/auGleriz} 260292 249 10nbT8ahi2 | 4ileti2ecias 
| i { { | i 
feo Ola52S be59eet5 on56fl2 7 $81{ 13 86 ieee ee oO! 
{ 1 { | | i 
feo 2 52 0 (eet OO ml bisa os eo hi 2a 95 3]) 30 or aot 
(he ed Le ed dhes ee be 
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Cards 9 through 36 would therefore be arranged as fcllows: 


{ CARD 9 contains 30 40 30; 


{ CAKD 10 contains 15 65 20; 


For { ° 
Approach { ° 
1 { ° 


{ CARD 14 contains 17 57 26; 

{ CARD 15 contains 41 36 23. 
The rest of the cards follcw the same arrangement as above; 
with CARDS 16-22 for Approach 2 CARDS 23-29 for Approach 3 
and CARDS 30-36 for Approach 4 in that order. 

In the above sequence of cards, CARDS 9, 16, 23 and 30 
contain the warm-up turn percentages from the Approaches 1, 
2, 3 and 4 respectively. 

If there are more sets of data of this form in the sense of 
period to be used and the percentages to be used; then 
different sets of the CARD 8 and CARDS 9-? are repeated, 
but again the number of such sets of data should not exceed 


the figure as represented in CARD 3. 


GoZzeo) Data, LYpe 3 

This is a combination of the first two types of data 
formats. In this type, estimates are made of percentages of 
the total volume cf vehicles that arrive at the various 
approach directions at the discrete time periods as may be 


specified, and at the same time estimate the turning 
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percentages of the estimated arrivals during the _ feriod 
specified. 

Here, it is possible to study the effect of periodic 
Changes both in arrival rates and turning probabilities by 
arbitrarily specifying arrival percentages (of the total 
hourly vclume) for discrete time intervals, and at the same 
time specifying the percentages of turning volumes for the 


particular period from all approach directions. 


Preparing input for Data Tyre 3 

The data format for Data Type 3 is shown in Tables C.5 and 
C.6. In order to manipulate data with the simulator as 
outlined above, the follcwing procedure of input data 


organization should be followed: 
(i) CARD 1 
should have a '3* in cclumn 4. 
(ii) to (vi) are the same as for Data Types 1 and 2. 
(vii) CARDS 9-? Correspond to CARDS 9-? Of Data Type 1. 


(viii) CARDS 16-? Correspond to CARDS 9-? Of Data Type 2. 
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Table C.5 
Sample Data Format 
Tctal Hourly Statistics 


im ~~ ai . a PEMA. Ae Ft = ee Se a 
{APPROACH j TOTAL { TURN VOLUMES | TURN FREQUENCIES | 
| VOLUME {| LEFT THRU RGT{ LFT THRU RGT | 
+-—--—_-_---++-—_--—___- -_-+---_-_-_--—___—__ —_ + -—--—-_ -— 4 
| 1 { 815 { 331 404 80, 0.40 0.49 0.09 | 
j { { { | 
{ 2 { 7188 ipt85af657 46], 0.01 0.83 0.05 | 
| { ( { { 
i 3 j 988 f 265 36680 43{ 0.26 0.68 0.04 | 
{ { { | { 
| 4 { 865 | 1S5a¢770 80] 0.01 0.89 0.09 | 
| ey he ee ae ee ee Dee J 
Table C.6 


Sample Data Format 3 
Pericdic Arrival and Turn Percentages 
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It is clear from the three types of data organization that 
Data Type 3 comprises Data Types 1 and 2 so that data types 
1 and 2 could be simulated using data type 3 by merely 
holding one of the generating characteristics (arrival 
volumes or turning volumes) constant. 

Of course, most parts of the input organization could 
be changed very easily by any one familiar with FORTRAN , by 
changing the appropriate format statementa. 

The user does net have to ccnform to the 6-column and the 4- 
column fields chosen for test data organization as described 


in the input organization section. 
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APPENDIX D 


ANALYSIS OF DIRECTIONAL FLOWS AT TRAFFIC CIRCLES 


This is a method developed by van Hotfen [50] for more 
accurate and precise analysis of traffic circle counts. 
Briefly, the method is as fcllows: 

Some selected movements are counted directly, and all 
other movements are calculated by solving simultaneous 
equations. Assuming that vehicles do not negotiate U-turns 
at straffiic circles; that is, vehicles entering the circle 
via the approach section of an artery do not exit from the 
circle via the exit section of the same artery; then 
treating cars and trucks alike (Vehicles), 12 equations can 
be set up involving 12 directional flows at any 4-arm 
traffic circle . 

The following abbreviations are made use of on the equations 
to follow: 


NE --- North-East traffic; 


NS --- North-South traffic; 
NW --- North-West traffic; 
WN --- West-North traffic; 
WE --- West-East traffic; 
WS --- West-South traffic; 
SW --- South-West traffic; 
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-- South-North traffic; 
-- South-East traffic: 
-- East-South traffic; 
-- East-West traffic; 


-- Fast-North traffic. 


Referring to Figure D.1 the indicated measured flows 


interpreted as follows: 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 


(7) 


(8) 


(9) 
(10) 
(11) 
(12) 


----> Vehicles entering circle from the North; 


<=> te ae aa ae '' West; 
=--—> ae ae ae ae 't South; 
ee ee ae ae ae ‘¢ East; 


---> Vehicles leaving circle at the North Exit; 

---> ae a qs ea t West “7 

os ae ae i] #¢ 06) 6 South os 

---> Vehicles which cross that point as indicated 
by the arrow in Figure D.1; 


---> Right-turning vehicles from the North; 


---> ee aa a9 't west; 
---> ae aa ‘# #8 South; 
> aa as it eet as it. 
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Figure D.! 
Measured flows to determine directional flows. 
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The 12 simultaneous equations that can be set up as a result 


of the above measurements are: 


(1)---> NE + NS + NW =a, ; 


(2)---> WN + WE + WS = az ; 

(3)---> SW + SN + SE = az ; 

(5) a= ote te N ag 3 

(5)---> WN + SN + EN = a, ; 

(6)---> SW + EW + NW Gye # 

(7)---> ES + NS + WS = a> ; 

(8)---> ES + NS + NE = ag ; 

(9)---> NW agus 

(10) --> ws aot 

(11) --> SE aes 

(t2)i=a> EN ajo 
Equations (9) to (12) are the only ones that are obtained 
directly from the field surveys. In addition tc the 
directly measured flows, the remaining 8 values are easily 
solved as follows: 


We set up 8 simultaneous linear equations with 8 unknowns in 


this manner. 


From 


From 


From 


From 


From 


(4) and ((9) ---~SMNEBtENS&= atzaa. S$; (A) 
(2) and (10) --> WN + WE = as-ajo ; (B) 
(3) and (11) --> SW + SN = Ber, ; (C) 
(Gard Geel eam | hott =a das daiar., (D) 
(S)eond (12) 5-—— oe WN +t oN) = as —ajo., (E) 
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Ws 


FLOMS(G)sanda (4 )t—--> SW + EWo= <a -—ac 3 ( F) 
From) manda (10)s==> ES t) NSe=i a7 ai: (G) 
From (8) -<----> ES + NS + NE = ag. (H) 
SOLUTIONS: 

Frome (G)7 and@(h)i---> NEV= vas (EStNS)) — 4, <a 4,,; 
HEB 90 A Ss a ae Aa Nua adie Qa ENTS 

Heatay lef (Eg) hohe Rat Ia dae POSE Oia NS 

BY ORE) Seige ce Be LINE Ee ti cine et Aree 

iGO Wee (Pt) coe Se Peo Ne eo eae TW; 

BO Mla (GC ) etn ne 2 SN >= a, -) diy = SN; 

ELON (E eos a — pln ee ZeNNG= ee Ape oN, 

BO OMGr (0) ee ZW Let Aone aioe a WN. 


Thus, knowing all the a's, the directional flows are easily 
computed from the selected measured flows. 

For example, the total flow from the four approaches 
and the correspcnding turning movements from the same four 


approaches are given as: 


Approach Volume Turn Movements 
From North 1 NE+NS+NW=a, NE NS NW; 
From West Z WNtWEtWS=an WN WE WS; 
From South 3 SW+SNtSE=az, SW SN SE; 
From East 4 EStEWtEN=a, ES EW EN. 


a) @5 a; a, are measured directly as part of the selected 
Counts tc bey nade: ala SW, SN, SE, ES, EW and so on, are 


obtained as solved above. 
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This method will fail however, when there is no 
guarantee for non U-turners at the circle, and it is 


applicable only te circles with 4 approaches and exits. 
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APPENDIX £ 


DATA COLLECTION AND ANALYSIS 


E.1 Data Collection for Peak-Hour Demand 

To be able to collect the data for the measured flows 
at a traffic circle inorder to determine the 12 directional 
flows of a 4-legged circle, 8 observers are needed for 
traffic counts. Fortunately enough, the model generates 
vehicle arrivals on the same basis as the counts made by the 
traLeic divisicn cf the City of Edmonton Engineering 
Department, that 1S, on approach basis; so that the Input 
Traffic Volumes used for the model generation routine were 
the peak-hour demand collected occassionally by the Traffic 
Division of the City of Edmonton Engineering and 
Transportation Department. 

The counts are one-hour traffic obtained for all 
approaches of a traffic circle. Traffic counts for each 
approach are stratified into left-turn, through, and fright- 
turn maneuvers. The data are for the morning and afternoon 
peak hours (739-830 AM and 439-530 PM respectively) for the 
following traffic circles in the City of Edmonton: 

1. 111 Avenue and Great Road; 

2. 114 Street and University Avenue; 

3. 125 Avenue and St. Albert Trail; 
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q. 114 Street and 72 Avenue. 


A brief description of each traffic circle is given in Table 


E.1 
Table E.1 
Summary of Study Circles 

— = as A. Ck 2 a oo ae 
{ TRAFFIC { TYPE | MAJOR { MINOR {| LOCATION | 
| CIRCLE { | { A { 
{ { | STREET { STREET { WITHIN | 
{ i | { { city { 
t-—__-—__-+--------—__--++ -----—-- —-+--_-—_-—_-++-—_--_-—_—_ 
{111 Ave | 4-Legged | { {Outlying | 
{ and | 2-way-2-way| Groat { 111 Ave |{ Shopping | 
{Groat Roadj 2 lanes in | Road { { Center | 
{ jall sections| | { | 
{ { { | { | 
{174 St. { 4-Legged { { {Outlying | 
{ and { 2-way-2-way{ 114 St. {Univ. AvejUniversity| 
{Univ. Ave | 2 lanes in | { { Area | 
{ fall sections| { { { 
| { | { { { 
{125 Ave {| 4-Legged { { {Outlying | 
{ and | Z-way-2-way|St. Albert{ 125 Ave |Shopping | 
{St. Albert{ 2 lanes in { Trail { { Center | 
{ Trail fall sections| | { | 
{ | { | { { 
{7114 St. { 4-Legged | { {Outlying | 
{ and | 2-way-2-way| 114 St. | =%72 Ave j[University| 
{72 Ave. { 2 lanes in | | { Area { 
{ jall sections| { { { 
ee Be A ae Se ee ee ee ee ee ee 


E.2 Data Collection to Establish Capacities 


of Circle Lane Sections 


The number of vehicles within a circle lane section is 
theoretically ccmputed from the specification of the length 
of a circle lane section and the effective vehicle length 


assumed in the systen. But due to the fact that the 


} HOTTAIOS | FOWIS ' 
: i 
i “iriw | TaRATS 
: 1772 | 
————-—— 
oatyieuo} I 
peiqgone | ova PTE 
iozas> | 


' 

' 

1 ! 

! iy 

eniyireo) 1 1 i 
Yiltewsvinviors .vintt .2e OTF 

no7vA =) ' 7 
: 
t 


' 


paky lang) 
pasqqede | svaé es 
tote | 


ytiedeviay) ava SY 
£924 


: 

! 
} 
' 
: 

: 
' 
‘ 
: 

i 
: 
1 
! 

: 
: 
: 
: 
: 
t 

L 


i} 
posylouo) : 
i 
' 


177 


vehicles in the circle lanes are always in motion (since 
they have absclute priority in the system), at is 
practically impossible to have all the vehicle spaces ina 
circle lane section simultaneously occupied. 

Data was therefore collected at some specific traffic 
circles to determine practically how many vehicles are 
within a lane section of a circle at any time during peak 
hour flows. One observer was adequate for the collection of 
this data, since attention was focussed on one circle lane 
section at a time first, attention was focussed on the inner 
circle lane section and for time intervals of 1 minute, the 
number of vehicles within the lane section was recorded. 
Thus after one hour, the maximum number of vehicles within a 
circle lane section during peak hour flow was established 


for both inner and outer circle lane sections. 


E.3 Data Collection for Gap Measurements 

The data collection for gap measurement at the circle 
was done in a simple less scphisticated manner. The circle 
lane sections to the left of a potential merger were divided 
into vehicle spaces for both inner and outer circle lane 
sections (see Figure E.1). The following procedure was 
adopted for the evaluation of gaps by inner and outer 


approach lane pctential mergers: 
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PYQure® Fa] 
Schematic of Circle depicting layouts for 
gap analysis studies. 
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For the evaluation of gaps by an inner approach 
vehicle, attention was fccussed on the first four vehicle 
positions of the inner and cuter circle lane sections to the 
left of the potential merger as shown in Figure E. Two 
observers were needed for this measurement. 

Whenever an inner approach vehicle accepted a gap to 
enter the circle, one observer would check the vehicle 
positions in the inner circle lane section which were 
occupied, and the other would do likewise for the outer 
circle lane section. This process was repeated for all 
vehicles entering the circle from the inner approach lane 


for a period of 30 minutes of the peak hour flow, and the 


results were tabulated as shown in Table E.2. 


Outer 

For the remaining 30 minutes of the peak hour flow 
‘attention was focussed on the first four vehicle positions 
of the outer circle lane section only. One observer was 
adequate for reccrding the evaluation of gaps by outer 
approach vehicles. Whenever an outer approach vehicle 
accepted a gap to enter the circle the occupied vehicle 
positions in the outer circle lane sections were checked 
only if an adjacent inner approach vehicle had not already 
accepted gap to enter the circle. That is, no frecord was 
kept for an outer approach vehicle taking advantage of gap 


analysis by an inner approach vehicle. 
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Table E.2 
Inner Lane Gap Analysis Summary Sheet 
I fe) 
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N---unoccupied vehicle space 
Y ---cccupied vehicle space 
Again, the process was carried out for all vehicles 
entering the circle from the outer approach lane for a 
period of 30 minutes of the peak hour flow, and the results 
were tabulated as shown in Table E.3. After analyzing the 
two tables for the peak hour data, it was deduced that the 
conflict zone for inner apfroach vehicles was the first two 
vehicle spaces of the inner and outer circle lane sections 
together with the last vehicle position of the potential 
merger's destination lane. Similarly, it was deduced that 
the conflict zone of the cuter approach lane was the first 
vehicle position of the pater circle lane section together 


with the last two vehicle positions of the potential 
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Table E.3 
Outer Lane Gap Analysis Summary Sheet 
OQ 
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N---unoccupied vehicle space 
Y---occupied vehicle space 
merger's destination lane. 

Such a gap analysis procedure is understandably not 
very accurate, but considering the amount and effort 
required in sophisticated field measurements, this crude 
method is deemed satisfactory for the accuracy required in 


the model. 
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